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We review some of the results obtained in the study of the production 
of electromagnetic radiation in relativistic nuclear collisions. We con- 
centrate on the emission of real photons and dileptons from the hot 
and dense strongly interacting phases of the reaction. We consider the 
contributions from the partonic sector, as well as those from the non- 
perturbative hadronic sector. We examine the current data, some of the 
predictions for future measurements, and comment on what has been 
learnt so far. 
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Electromagnetic radiation from relativistic nuclear collisions 3 

1. Introduction 

The study of matter under extreme conditions constitutes a rich field of 
intellectual pursuit and is a vibrant research area of physics. It is popu- 
larised to nonspecialists by indicating that such studies reveal physics that 
governed the early Universe (microseconds after the Big Bang) and also 
continue today to govern physics of compact astrophysical objects (neu- 
tron stars and black holes). But it is indeed more than that. Practitioners 
concern themselves with a variety of very specific and technically challeng- 
ing questions. When ordinary matter is heated to roughly a trillion degrees 
Kelvin, how does it respond? And what are the hallmark signatures of this 
response? What should one look for? When relatively cold matter is com- 
pressed to densities many times greater than that of normal nuclei, does it 
resemble something other than ordinary protons and neutrons? Quantum 
Chromodynamics (QCD) predicts under these extreme conditions, very far 
from the ground state, that matter will essentially change its properties to 
resemble a plasma of quarks and gluons (QGP)i^. What is the QCD phase 
diagram? 

Indeed, the understanding of QCD under extreme conditions of high 
temperature or large baryon density has progressed considerably in recent 
years. Confinement, in the pure glue version of QCD, is a property that 
can be associated with a definite symmetry whose status is probed by the 
value of the Polyakov loop (L). This symmetry is valid at low tempera- 
ture, but broken at high temperatures. In the limit of massless quarks, 
QCD is chirally symmetric, and that symmetry is valid at high tempera- 
tures and spontaneously broken at low temperatures. The order parameter 
there is the chiral condensate: {ipip). As a function of temperature, those 
order parameters are best studied on the lattice; although it is fair to say 
that the lattice has just started to venture into the finite baryon domain 
with any degree of quantitative assurance . An idea on the status of finite- 
temperature lattice QCD can be had by consulting Ref. 1, 2j. The order 
of the transition actually depends on the details of the parameters of the 
theory, as shown in the "Columbia plot" (Fig. but even in the absence 
of a robust prediction at finite temperature and baryon density, the lattice 
does provide tantalising clues of an eventually observable behaviour of the 
many-body nature of QCD. For example, it offers some support, together 
with effective models, to the discussions of a genuine tricritical point in the 
QCD phase diagra nP. 

The experimental realisation of hot and dense strongly interacting mat- 
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4 Electromagnetic radiation from relativistic nuclear collisions 




Fig. 1. The nature of the QCD phase transition as a function of quark masses, along 
with theoretical expectations from effective models. From 

ter in terrestrial accelerators has been accomplished through the study of 
relativistic nuclear collisions. Producing quark gluon plasma (QGP) is one 
of the premier goals of the Relativistic Heavy Ion Collider (RHIC) and 
its broader experimental program and represents already enormous chal- 
lenge; identifying the plasma and also studying its unique properties is yet 
something else, and has proven to be extremely challenging from both the 
theoretical and experimental points of view. Signatures of the QGP have 
been proposed: strangeness enhancement, suppression of the J/^p signal, 
effects of multiple collisions on the observed particle spectra, and electro- 
magnetic radiation are but a few. Unfortunately, while most of the proposed 
signatures are plausible and evidently do occur at some level, it has been 
difficult to refine the mostly heuristic arguments into really precise predic- 
tions. It is now clear that certainty will be attained in this field through 
the simultaneous analysis of complementary observables. This being said, 
probes that do not interact strongly have the definite advantage of suffering 
little or no final state interaction: they open a privileged window to the hot 
and dense phases of the reaction. The price to pay is a small production 
rate. 

We will discuss in this article the status of theory and the current picture 
relative to the experimental data, surrounding electromagnetic radiation as 
probes of strongly-interacting many-body dynamics. This work is organised 
as follows: a brief review of the formalism germane to the emission of real 
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Electromagnetic radiation from relativistic nuclear collisions 5 

and virtual photons from hot and dense systems is foUowed by a discussion 
of the radiation from the hadronic sector, then from the partonic sector. 
Some of predictions for current and future measurements are outhned, fol- 
lowed by a conclusion. The main thrust here pertains to models amenable 
to conventional experimental measurements: the high density low tempera- 
ture phase of QCD^will not be discussed, even though its electromagnetic 
emissivity has been investigated^. A goal here is to provide the reader with 
a snapshot of this rapidly evolving field by discussing some of the recent 
theoretical and experimental developments. In doing so, it is unfortunately 
impossible to do justice to the wholeness of the body of work in this excit- 
ing area: we shall be brief on some topics and refer instead to the literature 
and in particular to previous reviews. 

2. Radiation from Hadronic Matter 

The goal of this section of the text is to relate the spectrum of emitted 
radiation to some of the intrinsic properties of the strongly interacting 
matter. A thermalised medium is assumed, and the formalism below is 
developed in the one-photon approximation. 

Quite generally, the electromagnetic radiation from strongly interacting 
matter can be related to the imaginary part of the retarded in-medium 
photon self-energy, at finite temperature and densitjl^ . We sketch 

a derivatiorllil here, for the emission of real photons. Consider a transition 
i — /7, i.e. from some initial hadronic state i to some final hadronic state /, 
plus a photon of momentum = (cj, k) and polarisation e^. The transition 
rate is 

Rfr - (1) 

where t is the observation time, V is the volume of the system, and the 
^-matrix element is 



Sf,^{f\ j d^xJ^{x)A^{x)\i). (2) 
Jfi{x) is an electromagnetic current operator, and 



At'ix)^ -^={e'^-'' + er'^-'-'). (3) 

Using translation invariance, summing over polarisations, and using the 
integral representation of the delta function, one can write 
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6 Electromagnetic radiation from relativistic nuclear collisions 

x(/|J,,(0)|z)(*|J.(0)|/), (4) 

where the parts having to do with emission and absorption are evident. The 
net thermal rate is obtained by summing the above over final states and 
taking a thermal average of the initial configurations. 

It is useful at th is point to define some finite-temperature current- 
current correlator!^ 

i 

•' i 

f^Uk) - / d4xe*^5^(z| [A,ix),AAO)] K> e-'^^'/Z . (5) 
i 

The last correlation function is a retarded correlation function. The above 
all involve the electromagnetic current operator in the Hcisenberg picture 
and arc written in the grand canonical ensemble, where Z is the partition 
function. Assuming translational invariance, the first two can be rewritten 
together as 

= ^(2^)4<5^p/ -p, ±p)(*|A^(0)|/) (f\A,iOMe-^^'/Z. (6) 

Clearly, is involved with absorption of radiation, whereas deals with 
emission. Only, the latter case is treated here. Further defining a spectral 
density p^^, = f>^ ~ /<^, one may first show that 

and also that 

With the above elements, one can finally write — 2Im/^^. The emis- 
sion probability is related to the imaginary part of the finite-temperature 
retarded current-current correlation function. In the one-photon approxi- 
mation (i.e. to lowest order in e^), the time-order ed c urrent correlator is 
the one-particle irreducible photon self-energy, Il^Ji^. Putting all of this 
together, the differential rate for emitting real photons is 
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Electromagnetic radiation from relativistic nuclear collisions 7 

The proof is easily generalised to the case of lepton pair emission (the lepton 
mass has been set to zero): 

' (10) 



Thus, the electromagnetic signal emitted during a nuclear reactions can 
be related to a quantity that is linked to properties of the medium itself. 
As it shall be seen, the retarded self-energy (or alternatively the current- 
current correlator) will be modified in a strongly interacting environment. 
Also, the current-current correlator is calculable only perturbatively, unless 
a specific model is available: this again testifies to the importance and the 
value of electromagnetic measurement in nuclear collisions. They indeed 
open a window to the hot and dense phases of the reaction, and those 
regions can't be probed directly by other means. 

2.1. The Low Dilepton Invariant Mass Sector 
2.1.1. A Baseline Calculation 

As a calculation to set the scale of the physical processes under consid- 
eration, it is useful to consider first the following question: what is the 
magnitude of the radia tion emitted by a hot gas of mesons? Specialising 
to the lepton pair sectoJ^l this problem is briefly summarised here. Using 
relativistic kinetic theory, the lepton pair emission rates can be calculated 
with the help of effective interaction Lagrangians. The parameters of those 
effective Lagrangians are fitted to radiative decays measurements using vec- 
tor meson dominance (VMD). Specifically, the "calibration" reactions are: 
p 7r7, ^ K^-i, K *o (if*°) -> /s:0(X*°)7, uj -> 7r"7, p° r]j, 

r( — > p°7, r( — > 017, — > 777, — > 77'7, and 7r°7. Note that the usage of 
relativistic kinetic theory here is tantamount t o ev aluating the finite tem- 
perature photon self-energy at the one-loop level^. The rate for ah — > e^e^ 
is given by 

p _ ^ \r f '^^P" '^^P'' f f 

rtab^e+e- J 2Ea{2nf 2Ei.{2^f 2E+{2TTf 2E^{2^f ^"^^ 

x\M\'{2^f5\pa+Pb-p+-P-), (11) 

where the /'s are appropriate distribution functions, and M is an overall 
degeneracy factor dependent upon the specific channel. The PP e+e^, 
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8 Electromagnetic radiation from relativistic nuclear collisions 

PV — > e"*"e~, and VV e^e~ reactions can be obtained from the radiative 
decay ones through crossing symmetry. The sum of them is shown in Fig. [51 
Also shown is the contribution from the tt+tt^ e+e^ reaction. This 
channel has often been considered as the sole source of hadronic dileptons 
in early calculations, owing mainly to multiplicity arguments. Even in this 
incoherent sum approach, one can see that this assumption is badly violated 
in the low mass region. 



10^ 



> 



a 10"' 



X) 




IO- 
CS 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

M (GeV) 

Fig. 2. The rates for dielectron emission from an incoherent sum of meson reactionJi^ 
The three sets of full curves are the net rates at a temperature of 100, 150, and 200 MeV, 
from bottom to top. The dashed curves represent the tt+tt" — > e^e~ contribution only. 



2.1.2. Spectral Density Calculations 

The full power of the formalism derived in section [3 reveals itself when 
it is made clear that the electromagnetic radiation rate is related to the 
in-medium vector meson spectral density: this quantity is of course not 
measurable directly. The spectral density is related to the imaginary part 
of the full propagator, by definition. It appears when the electromagnetic 
current operator in the retard ed self-energy is expressed as the field op- 
erator, through We illustrate the spectral density approach 
here, by describing a calculation again done in the baryonless regim^2D 
One may start with the model for the p-meson in free space employed 
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Electromagnetic radiation from relativistic nuclear collisions 9 

previously in Refs. |15[ 1161 [T7 ] . Based on the standard pinT interaction 
vertex (isospin structure suppressed), 

Cp-„T, = 5p^7r TT P^T^Pn , (12) 

(p^ : pion momentum) the bare p- meson of mass is renormahsed through 
the two-pion loop including a once-subtracted dispersion relation, giving 
rise to the vacuum self-energy 

t%^{M) = j ^ ^pUvf GlJM,p) , (13) 
with the vacuum two-pion propagator 

^-(^^'^)-^ M^-(2.!(p))^+., ^ -^iP) = V^<^, (14) 
and vertex functions 

VpTT^ip) y'^ gp^TT 2p Fp^^{p) , (15) 

involving a hadronic (dipole) form factor Fp^^ 1121 Resumming the two-pion 
loops in a Dyson equation gives the free p propagator 

D°p{M) = [Af2 - (^0)2 _ sO_(Af)]-i , (16) 

which agrees well with the measured p -wave tttt phase shifts and the pion 
electromagnetic form factor obtained within VMD. 

To calculate medium corrections to the p self-energy in a hot meson 
gas, one can assume that the interactions are dominated by s-channel res- 
onance formation. It is then possible to use the formal relationship that 
relates the self-energy to the forward scattering amplitude, integrated over 
phase space at finite temperature. A field-theoretic derivation of this con- 
nection can be found in |22| . At moderate temperatures relevant for the 
hadronic gas phase, the light pseudoscalar Goldstone bosons P = tt,K are 
the most abundant species. The various resonances in pP collisions can 
be grouped into two major categories, namely vector mesons V and axial- 
vector mesons A. The effective Lagrangians that regulate the interactions 
among all those species have their parameters chosen such that the mea- 
sured hadronic phenomenology are reproduced'*''^. This statement also holds 
true for hadronic form factors. In order to calculate the p spectral density 
at moderate temperatures, one includes the hadronic fields appearing in 
Table ^ The interaction vertices being completely determined, one may 
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Fig. 3. The real and imaginary parts of the polarisation-averaged p self-energy (lower 
and upper panel, resnectively ) . The different channels are labelled explicitly and are 
explained in the texi=JJ Note that the tttt channel is absent for the sake of viewing 
clarity. 



Table 1. Mesonic resonances R with masses m/j < 1300 MeV and substantial 
branching ratios into final states involving direct p's (hadronic) or p-like photons 
(radiative). Taken from Ref. 1211 . 



R 




rtot [MeV] 


ph Decay 


V\ [MeV] 


[MeV] 


w(782) 


0-1- 


8.43 


pTT 


~ 5 


0.72 


/ti(1170) 


0-1+ 


~ 360 


pTT 


seen 


7 


ai(1260) 


1-1+ 


~ 400 


pTT 


dominant 


0.64 


7^1(1270) 


il+ 

2 ^ 


~ 90 


pK 


~ 60 


7 


/i(1285) 


0+1+ 


25 


PP 


<8 


1.65 


7r'(1300) 


1^0- 


~ 400 


p'K 


seen 


7 



first calculate the in-medium p self-energy, th en t he complete propagator 
in terms of its longitudinal and transverse partW 

The real and imaginary 

parts of the in-medium p self-energy are shown in Fig. |3 Each curve is 
labelled according to that species which interacts with the p. It is instruc- 
tive to observe that the imaginary parts all add, while there is a significant 
amount of cancellation of the real parts. The first effect creates a sizeable 
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Fig. 4. Imaginary part of the p propagator in the vacuum (solid curve) and in a thermal 
gas including the full in-medium self-energies for fixed three-momentum q = 0.3 GeV at 
temperatures T = 120 MeV (long-dashed curve), T = 150 MeV (dashed curve), and T 
= 180 MeV (dotted curve). The figure comes from Ref. |21| . 



width for the in-medium p, while the second determines the in-medium 
mass, which appears to be only slightly modified. Those different aspects 
are again seen in the representation of the imaginary part of the in-medium 
p propagator, shown in Fig. 0J 

The calculations of the in-medium vector meson spectral densities 
clearly show the richness of the many-body problem under scrutiny. The 
power of the techniques described above becomes evident when they are 
combined with dynamical models and confronted with experimental data. 
This story is well chronicled in j^SI 1^ , and in references therein. See also 
|25j . The current situati on can be summarised by writing that the low dilep- 
ton invariant mass 

dataPEZI 

can be understood in terms of in-medium 
modifications of vector spectral densi ties. Thes e data can not empirically 
exclude, however, other inter pret ationJ ^ I "^'^ ! "^'^ ! This unfortunate situation 
still prevails at a lower energy"^ as shown in Fig.El The sources there are: 
free hadron decays without p decay (thin solid line), calculation with a vac- 
uum p spectral density (thick dashed line), dropping in- medium p inas m 
(dash-dotted line) , and with a medium- modified p spectral density^ (thick 
solid line) . Note that but the persistence of the dilepton excess at lower en- 
ergies does support a baryon density-driven effe ct. Other suggestions to 
resolve the different models require high statisticJ^ the final analysis of 
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0.2 0.4 0.6 0.8 1 1.2 



m^, (GeV/c') 

Fig. 5. Dilepton spectrum from Pb-Au collisions at 40 AGeV/c. See the text for the 
meaning of the different curves. Data are from 1^. 

the CERES 2000 data with the TPC is eagerly anticipated. 

2.2. The Dilepton Intermediate Invariant Mass Sector 

Intermediate mass lepton pairs have traditionaUy been a focus of inter- 
est, as their spectrum has been suggested early on as a signature of the 
quark-gluon plasma . In relativistic nuclear collisions, measurements have 
been carried out at SPS energies by the IIELIOS-3 and the NA38/NA50 
collaborations in the lepton pair invariant mass range < M < mj/^. 
Both experimental collaborations have observed significant enhancement 
of dilepton yield in this region for central S + W and S + U colli- 
sions as compared to those i n prot on-induced reactions (normalised to the 
charged-particle multiplicity)"^'""'. Chronologically, HELIOS-3 reported on 
the intermediate-mass enhancement first. This experiment was designed to 
study virtual photons in the dimuon sector at low transverse mass. In this 
way, dimuon production was studied from threshold up the J /ip mass over a 
wide range inpx- A good summary of this experimental situation is shown in 
Fig. El Several explanations of the intermediate mass dimuon enhance ment 
have been put forward. Those include additional production of cc pair 
secondary Drell-Yan emissiort^^, and charmed meson rescattering^"! Note 
that in principle, all of those effects can coexist. However, such a global 
modelling has not been done, and we thus discuss them separately. It is 
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M (GeV) 

Fig. 6. Comparison of lepton-pair yield divided by the multiplicity of charged particles, 
in p + W and S + W collisions at 200 A GeV/c. The data are from l27l . 



fair to say that the first of those mentioned above is still admitted by cur- 
rent experimental data ' , even though some calculations are slightly less 
optimistic . The role of the last ingredient on the list has not been found 
to be large enough to account for the experimentally observed excesJ^^^^. 

Another class of approaches consists of quantitative evaluation of ther- 
mal dilepton sources. Those may be from the hadronic, confined sector of 
QCD, and/or from the quark gluon plasma itself. One such model is de- 
scribed below. Recall that thermal hadronic sources have been shown to be 
crucial in the low mass sector. It is therefore legitimate to ask how high 
in invariant mass is the extent of the virtual photon radiation from those 
sources? Those concerns are carried to their logical conclusion in what fol- 
lows. 

In the intermediate invariant mass region, relativistic theory estimates 
indicate that the following microscopic channels are relevant: tttt — > , 
TTp^l+r, Tiuj l+l', TTfli ^ KK -^l+l-, and KK* + c.c. l+l' 

|4()j . Apart from sheer coupling constant values, the importance of those 
contributions stems simply from considering energy scales involved and 
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14 Electromagnetic radiation from relativistic nuclear collisions 

from phase space arguments. This combination of couphng constants and 
phase space is effective in m aximising a particular contribution from the 
TTfli channel, for examplJ^. 

Calculations of dilepton-emitting processes in the intermediate invariant 
mass region follow similar steps to those in the low mass sector. Effective 
Lagrangians are used, together with VMD, and the coupling constants and 
possible form factors are fitted to measured strong decays and electromag- 
netic radiative decays. Only, in the intermediate mass domain an extrap- 
olation is required. The strong decay widths set a scale that is typically 
an order of magnitude below the mass region of interest: 1 GeV < M < 
3 GeV. The radiative decays are even smaller, owing to the size of Ofcm 
|42| . The required extrapolation is then vulnerable to off-shell effects. Put 
another way, there is a risk of uncontrolled growth of form factors since the 
application region is far removed from the region where the empirical fitting 
was realized. Indeed, different Lagrangians known to agree in the low mass 
sector generating dileptons were found to differ significantly in their predic- 
tions of emission rates for intermediate mass lepton paira^. Fortunately, 
there exists a wealth of data for e~^e~ —* hadrons, exactly in the invariant 
mass window relevant for this application. Those data can thus be used to 
extract an effective form factor for the inverse reactions. Alternatively, they 
may also be used to extract spectral densities: this point will be discussed 
later. As an example, consider e'^e^ tt+tt^, which has been measured 
with high accuracy evidenced by the data shown in Fig. [3 The cross 

section for this reaction can be written as 

^(e+e--.7r+7r-) = ^^|F.(Af)p, (17) 

where k is the three-momentum in the two-body rest frame, M is the lepton- 
pair invariant mass, and F^^ is the time-like pion electromagnetic form fac- 
tor. With these data, one can extract F^^, and then use it in the calculation 
of the dilepton-producing reaction tt+tt^ —>■ l^l^: 

a(7r+vr- ^ /+r) = — _|F.(M)p 1 _ h + ^ , (ig) 



where mi is the lepton mass. A similar procedure can be followed for other 
channels. Another example appears in Fig. |S1 where the time-like elec- 
tromagnetic form factors for the kaon systems have been extracted from 
electron-positron annihilation data. Both processes introduced above are 
of the pseudoscalar-pseudoscalar type. For the pseudoscalar- vector class, in 
the invariant mass region of interest, vrp l^l^ , KK* + c.c. l^l^ , and 



February 1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3v8 



Electromagnetic radiation from relativistic nuclear collisions 



15 



lO" 



1000 



^ 100 



10 



0.1 




0.5 



Fig. 7. The cross section for e+e — > Tr+Tr . The solid curve is base d oil the model 
of Ref. 1461 . Th e experimental data are from the OLYA coIIaboratioJSl the CMD 
coUaboratior^ and the DM2 collaboration EHl 



■KLo — > /+/~ are included. The first two processes effectively involve three 
pions, while the third one involves four pions. Note that in a transport 
approach, a process involving three or more pions in the initial state can 
only be described as a two-step process with an intermediate resonance. 
The first two channels above have been studied in Ref. [SII. The effective 
form factors one extracts are shown in (41^. Details about the ttuj channel 
are gotten from the study of four-pion final states. Using a Wess-Zumino 
VVP interaction Lagrangian, one finds 



9M 



(19) 



in the limit of vanishing lepton mass. The form factor may be parametrised 
in terms of three isovector p-like vector mesons, p(770), p(1450), and 
p(1700jSIl 



Ft:uj{M) - V.gyTrwT— 2 T7o\ ■■ 

[rriy — M^j — imyl y 



(20) 



V 
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Fig. 8. The kaon electromagnetic form factor. The solid and dotted curve are for the 
charged and neutral kao ns. r espectivelv. The svmbols are for the charged kaon data from 
the CMD-2 collaboratioiE2| the DM2 coUaboratioiEU and the OLYA coIlaboratiorEHl 



This form is then used to fit the experimental data from the ND and AR- 
GUS collaborations. The result is shown in Fig. |21 

In the pseudoscalar axial-vector channel, we shall consider mainly 
TTfli l^l^ , which is in effect a four-pion process. Considering here again 
the reaction where the lepton pair constitutes the initial state and the 
hadrons the final, one can attempt an extraction of an effective form factor. 
Some previous thermal rate calculations indicate that this specific channel 
is particularly important in the intermediate mass regioJ^, even though it 
is difficult to calculate reliably a specific signal using effective Lagrangians. 
This fact owes mainly to off-shell effectJ^. One can pick a model that yields 
adequate hadronic phenomenology on-shcU, and then extrapolate to the in- 
termediate mass sector with the help of experimental data. Using a chiral 
Lagrangian where the vector mesons are introduced as massive Yang-Mills 
fielda^ one may derive the following cross section 



IV 
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2ml 



(21) 



1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3v8 



Electromagnetic radiation from relativistic nuclear collisions 17 




M (GeV) 



Fig. 9. The cross section for e+e — > tt^oj. The so lid c urve is described in the text. The 
experimental data are from the NrEHand ARGU^^ collaborations. 



where 7i is a nontrivial function of coupling constants, masses and mo- 
menta, is the magnitude of the pion momentum in the centre-of-mass. 
The issue of the electromagnetic form factor \F.^ai can be settled, at least 
in principle, by analysing e"'"e~ — > 7r+7r~7r+7r~ and e^e~ — > 7r+7r~7r'^7r'^ 
data. Although many such analyses have been carried out, an unambigu- 
ous result is still elusive, as many other intermediate states may contribute 
to the same four-pion final state. Several scenarios have been considered, 
and a discussion appears in '40' . What is probably a conservative estimate 
is highlighted here. The DM2 collaboration has determined the cross sec- 
tion i^e+e- ^nai using a partial wave analysis One may extract 
an effective form factor from these data, see Fig. and then carry out an 
analysis for crural -tifi using detailed balance. 

Even in a careful analysis of the relevant intermediate invariant mass 
dilepton reactions, some concerns remain. These mainly stem from the need 
to account for all the sources of electromagnetic radiation. In kinetic the- 
ory approaches, risks exist of double-counting and possible omissions. In 
an attempt to bypass those, an approach which allows for a nonperturba- 
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Fig. 10. The cross section for e+e~ — > irai. The open circles are the experimental data 
from the DM2 collaboration using a partial wave analysiiSl The solid curve is a fit to 
the data. 



tive treatment of the strong; interaction and avoids a detailed enumeration 
of reactions was deviscd'^'^. The dilcpton emission rate is interpreted in 
terms of spectral functions of liadronic currents, tabulated from low energy 
e~^e~ annihilation reactions and from t lepton decays. A differential rate 
expression, obtained in the chiral (m^ — > 0) limit, readJ^ 



dR 



4a2 

= MTKi{M/T) 

2tt 



{p''{M)-p^{M)) 



(22) 



where T is the temperature, e — T^/6F^, M is the dilepton invariant mass, 
and the superscripts on p denote the electromagnetic, vector, and axial 
spectral functions, respectively. These spectral distributions are displayed 
in Fig. 111! Using the spectral functions to generate the lepton pair emis- 
sion rate, a comparison with the rates obtained via a summation of mesonic 
reaction channels is shown in Fig. 1121 To summarise, the contributing chan- 
nels producing lepton pairs in the invariant mass range 1 GcV < Af < 3 
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Fig. 11. The spectral functions p'^"'{s), p^{s), and p^{s), as compiled from 
nihilation and t decay data^i^ 



GeV have been found to correspond to the initial states tttt, irp, ttlu, rjp, 
pp, TTQi, KK, KK* + c.c. [S^. The detailed channel-by-channel assessment 
clearly accounts for the net signal yielded by the "global" spectral function 
analysis. 

In order to compare with experimental data, the rates must be time- 
integrated in a model that is also compatible with other measured ob- 
servables, hadronic or otherwise. Furthermore, a precise simulation of the 
detector acceptance and resolution is necessary. An approach that incorpo- 
rates both aspects is described presently. A class of models that produce 
time-evolution scenarios is that of hydrodynamic models. Specifically, the 
assumption is that, at SPS energies, a plasma is produced at proper time 
tq. Assuming isentropic expansion, the temperature and pr ope r formation 
time can be related to the measured differential multiplicity^ 

dN/dy is the measured particle rapidity density and a = 42.257r^/90 for a 
plasma of massless u, d, s, g partons. Once the transverse area At is known 
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Fig. 12. Net dilepton production from a gas of mesons at a temperature of T = 150 
MeV, as a function of dilepton invariant mass. The solid curve is the sum of the hadronic 
channels discussed in the text and in the references. The data points are from 157) . 

along with dN/dy, the above relation links To with tq. Enumeration of the 
model premises is completed by the statement that the plasma is assumed 
to undergo a boost-invariant longitudinal expansion and an azimuthally- 
symmetric radial expansion, with a transition to a hot hadronic gas con- 
sisting of all hadrons having M < 2.5 GeV, in thermal and chemical equi- 
librium at temperature Tc- This makes for a rich equation of state. Once all 
parton matter is converted into hadronic matter, expansion continues un- 
til a kinetic freeze-out temperature Tp is reached. Those steps are generic 
in hydrodynamic calculations. Note that during the evolution, the speed 
of sound in matter is consistently calculated at every temperature that is 
input int o t he equation of state and needed to solve the hydrodynamic 
equationJSSI Additional details about setting up the initial conditions for 
the hydrodynamic evolution can be found in Ref. The same reference 
also shows the result of hadronic spectra calculations with the hydrody- 
namic approach. 

It is vital to account for the finite acceptance of the detectors and for 
their resolution when comparing the results of theoretical calculations with 
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measured experimental data. In the case at hand, those effects are indeed 
important in the NA50 experimenJ^. One approach to this problem in the 
past has been to model approximately and analytically the acceptancJ^^^^. 
While this can be readily implemented, a legitimate doubt can subsist about 
the accuracy of the experimental representation, especially in regions where 
edge effects might be important. In order to circumvent this problem, a nu- 
merical subroutine developed to reproduce the NA50 acceptance cuts and 
finite resolution effects in the measurement of muon pairs in Pb + Pb colli- 
sions at the CERN SPS was uscq*^^. Thus, the invariant mass distribution 
of lepton pairs is computed in the hydrodynamic model, and then the pairs 
are run though the numerical detector simulation. The normalisation is de- 
termined by a fit to the Drell-Yan data using the MRSA parton distribution 
functions, as in the NA50 analysis. In order to get a pr distribution, the 
dN/dA'P estimates for Drell-Yan were supplemented with a Gaussian dis- 
tribution in pi*^, this very closely reproduces estimates obtained by the 
NA50 collaboration. The resulting invariant mass distribution is shown in 
Fig. El The pr distribution is also computed. It is shown in Fig. In 




2 3 4 5 

M (GeV) 



Fig. 13. The calculated dimuon invariant mass distribution, after correctin g fo r the 
detector acceptance and resolution. The data are from the NA50 coUaboratioiiSsl The 
Drell-Yan and thermal contributions are shown separately, as well as those coming from 
correlated charm decay and from the direct decays of the J/i/i and i/"'- 



both cases, good agreement with the experimental data is clearly achieved. 
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Fig. 14. The dimuon transverse momentum spectrum, after aceounting for the detector 
effects. The data and the different curves are from the same sources as in Fig. 1131 

As this point it is appropriate to consider the foUowing question: which 
initial temperature is demanded by the intermediate invariant mass dilepton 
data? A critical and quantitative assessment of this issue can be obtained 
by examining a linear plot of the lepton pair mass spectrum in the region 
under scrutiny. This is shown in Fig. 1151 From this figure it is clear that 
the best fit is provided by tq = 0.2 fm/c, and that the second best (less 
than two standard deviations away for most of the data points) belongs 
to To — 0.4 fm/c. In terms of initial temperatures, those correspond to 
To « 330 and 265 MeV respectively. A conservative and reasonable point 
of view is that it is probably not fair in such a challenging and complex 
environment as that of ultrarelativistic heavy ion collisions to ask for an 
agreement that is better than two standard deviations, considering all of 
the inherent uncertainties. The quark matter contribution (as modeled by 
qq annihilation) is w 23% for tq = 0.2 fm/c, and ~ 19% for tq = 0.4 fm/c, 
around a lepton pair invariant mass of 1.5 GeV. 

Focus so far has been placed on high multiplicity data only. However, 
to extend the hydrodynamic model to non-central events and to properly 
treat the azimuthal anisotropy is not a simple task. However, one can get an 
approximate estimate of the centrality dependence by ignoring the broken 
azimuthal symmetry and by approximating the region of nuclear overlap 
by a circle of rad ius R « 1.2 (iVpart/2)^/^, where iVpart is the number of 
participantJS^ESl. A centrality-dependence is generated thusly and shown 
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Fig. 15. A linear plot of the net dilepton spectrum in the intermediate mass region. 
The three solid curves correspond to formation time to= 0.2, 0.4, 0.6 fm/c, from top 
to bottom, respectively. The data are from i6fl . The thermal contribution and that for 
hard processes are shown separately. 



in Fig.^l It is seen that the agreement with the measured data is quite good 




30 50 70 100 200 300 



Np„t 

Fig. 16. Centrality dependence. The data represents the "extra charm" yield (as char- 
acterised by the NA50 collaboratiorl^3I) needed to describe the intermediate mass dimuon 
data. The solid curves are from the sources discussed in the text. 




and that this approach gives a fair description of the centrahty dependence 
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of the excess dilepton measurement. 

What has been learnt from this exercise is that many-body channels, 
important in the low mass region, also still play a vital role in the interme- 
diate mass domain. This unity is satisfying. It also happens that the global 
dynami cal behaviour of the electromagnetic radiation can be empirically 
modeled^. Perhaps more importantly, the case where the data shown in 
this subsection are interpreted as the signature of a charm excess no longer 
appears to be very compelling. The findings described here are in agreement 
with previous calculations of dilepton radiation at this mass scal e ' 'I As 
mentioned previously, the portion of the signal that emanates from the de- 
confined sector is around 20%, a figure that is unfortunately too small to 
support convincing claims of a QGP presence, once all uncertainties are 
factored in. 

2.3. Photons 

2.3.1. General Strategy 

Real photons differ from previously discussed dileptons in a couple of im- 
portant ways. First, they are on-shell and thus cannot be accommodated 
kinematically with two-hadron annihilation processes. From the beginning 
then, the only one-loop contributions to the current-current correlator (or 
retarded self-energy) are the hadronic radiative decays tt*^ — > 77, 77 ^ 77, 
and uj — > TT^^. Although the oj lifetime is ^ 23 fm/c and so that channel 
could be considered a thermal source (just barely), the others are clearly 
non-thermal sources with lifetimes much longer than the fireball. Such con- 
tributions must be considered in the overall yield, but are outside of the 
scope of the present discussion focusing on thermal emission. The one- 
loop contributions therefore play a smaller role in photon production as 
compared with dilepton production and consequently the discussion begins 
seriously here at two loops. Second, and something of a technical point, 
is that real photons have only two polarisation states over which to sum, 
rather than three as was the case for virtual photon propagation. In the 
same spirit then as the dilepton case, real photon emission from resonance 
hadronic matter can be most systematically studied by beginning with the 
photon self energy at two-loop order and working upward. The occupation 
numbers for the internal lines pay always a Boltzmann penalty and it is 
therefore natural to begin with the lightest species and with the minimum 
number of hadrons. As a aside remark, since the evaluation of two-loop 
topologies at finite temperature is technically challenging (although not 
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impossible), many photon rate calculations rely on a kinetic theory ap- 
proach. 

Since the energy regime is both relativistic and nonperturbative in terms 
of QCD degrees of freedom, it is commonplace to use effective theories for 
the composite hadron dynamics. Typically one starts with an effective La- 
grangian with a large enough flavour symmetry to account for the lightest 
and relevant species. As a general rule, the pions are most important, fol- 
lowed by rho and so on, simply owing to increasing mass. Quantum numbers 
also play a role in terms of spin states and isospin states governing densi- 
ties, and so one must be systematic. With interactions under some control 
relative to chiral symmetries, gauge invariance, conservation requirements 
of various sorts, one uses cutting rules on the two-loop self-energy diagrams 
in order to generate a list of reactions of the type ha + hb ^ hi + j and 
ha ^ hi + h2 + 7. 

2.3.2. Establishing the Rates 

The above mentioned strategy was first taken by Kapusta, Lichard and 
Seiberl^i where ir-p and light meson dynamics were investigated. The dy- 
namics were modeled with 

C = ^\D,<S>{' - ml |$p - ^p.^p^'^ + \mlp,,p^ - -F^^F^^ . (24) 

Coupling of the rho and the photon to pions was accomplished with the 
covariant derivative Z)^ = — ieA^^ — igpPfj,. The charged and neutral 
pions are embodied in the complex pseudoscalar field the vector rho 
and photon field strength tensors are respectively pp,i, = df^pu — duPn and 
F^v = d^Aiy — d,jAf^. Calibration is done by fitting the p tt+tt^ decay 
rate with the choice gp = 2.9. 

The specific channels studied in Ref. [fifij were dubbed annihilation 
TT^TT" — !■ p^ and "Compton scattering" irp —f ttj, and finally, neutral rho 
decay p — s- 7r+7r^7 (essentiall y t he finite temperature analog of the vac- 
uum process studied by Singep^^. Since the 77 meson mass is intermediate 
between pion and rho, its Boltzmann penalty is less than rho's. Its effects 
were also considered by including the channels 7r"'"7r~ — > 777, tt^ — > n^j. 
Owing mostly to coupling strengths (or weaknesses), these channels were 
found to be less important as compared to the purely tt and p channels by 
more than an order of magnitude. Finally, in this initial study of photon 
production, the channel tt+tt" — > 77 was included, though it was seen to 
contribute very little. 
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The matrix elements for all the processes enumerated above are included 
in Ref. and will therefore not be repeated here. The energy dependent 
invariant rate for producing photons is then obtained by folding in Bose- 
Einstein (enhanced, if final) hadron distribution functions and Lorentz in- 
variant phase space. For instance, for the channels Pa + Pb ^ Pi + P-y , one 
has 

xr /" I A-^ |2/o ^4r4/_ , „ „ „ \ d^Pa 



{2Tr)^2Eb ■"' {2TTf2Ei ' ' " (27r)32' 
where A/" is the appropriate degeneracy factor counting the states. 

The resulting rates have been established numerically. However, analyt- 
ical parametrisations valid for 100 MeV < T < 200 MeV and 0.2 GeV < 

< 3 GeV have been proposeo ^. 



dR 

E^-jj—inn p7) = O.OTITT^-^^^ exp(-0.7315/r + 1.45/^/-^.^ - E^/T), 
Ey-^i^P ^ ^l) = T'-' eM-y{2TE^f/^ - E^/n 
E^^{p ^nn^) = 0.1105 T4-283iJ-3-076+0.0777/T exp(-1.18i?^/T).(26) 

In these expressions, T is the temperature, E^ is the photon energy; both 
must be reported in GeV. The numerical constants have appropriate units 
in each case and the numerical values out in front have units fm~^GeV~^. 
The numerical results are shown in Fig. 1171 

The rate for the u! 7r*'7 is also included in Fig. El It can be written 

as'' 

p r , 3m^r(a;->^"7) 7 ,„ , . 

X [1 + fBE {E^ - E^)] (27) 

where -Emin — ''n^{E'^ + Eq)/2Ej Eq, and Eq is the photon energy in the 
rest frame of the uj meson. 



^The process p — ► tttt^ was slightly miscalculated in Ref. 1661 owing to an omission 
of a Lorentz-boost factor. The parametrisation of the process published in Ref. 1681 
is therefore not optimal. A slightly different parametrisation is proposed here which 
correctly accounts for covariance effects. 

''Note that Eq. (54) of Ref. 1661 has an incorrect Lorentz-boost factor. The formula has 
been corrected here and reported in Eq. 1271 . 
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Fig. 17. Photon emission rates from Ref. |681 plus an estimate including the ai coher- 
ently. The u! 7r7 rate results from Eq. 1271 . Temperatures are fixed at T = 150 MeV 
and T = 200 MeV. 



2.3.3. Refinements 

Soon after these initial rate calculations were done, Xiong, Shuryak and 
BrowiH^ pointed out that the ai meson would have an important effect 
on the Tip 7r7 channel. It is a resonance in the irp sector with pole 
mass roughly matching the average y/s in the fireball and with a rather 
large width ^ 400 MeV. An interaction for the aiirp vertex was proposed, 
stemming not from some symmetry argument but rather, a vertex hav- 
ing minimal momentum dependence and still respecting gauge invariance. 
The idea was to set the scale as simply as possible for the contribution 
from this process. From the strong-interaction vertex, vector dominance 
was employed to subsequently describe radiative decay. The ai was indeed 
found to be important when studied this way — even dominating the other 
exchanges. 

However, by itself the s-channel oi diagram does not carry complete 
information on the overall strength of the Compton process. Instead, a co- 
herent sum of pion exchange and ai exchange was necessitated. SonJ^ 



February 1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3v8 



28 Electromagnetic radiation from relativistic nuclear collisions 

later carried out a study of photon production using a Chiral Lagrangian 
with vector and axial-vector fields introduced as dynamical gauge bosons 
of a hidden local symmetry. The role of the ai once again seemed to be 
important. However, Song found that fitting the relevant masses and cou- 
pling strengths in the model allowed two solutions, i.e. two parameter sets, 
and was therefore not able to uniquely identify an emission rate from this 
channel. It was later found, while studying dilepton emission with the same 
model, that parameter set II was the more r elia ble one in terms of its ability 
to match the observable hadronic quantitiel^. In particular, the D/S ratio 
in the scattering amplitude for oi np was more closely respected with 
parameter set II. One thus observes a diminished presence of the oi meson 
(as opposed to estimates cited above) in the final rate. A modest enhance- 
ment of 20%, attributable to the pseudo- vector, is illustrated in Fig. El An 
update on this is forthcoming^. A study with the hidden local symmetry 
approach also points to a reduced role of the oi mesorn^. 

Strange particles have been found to be only marginally important in 
the literature up to now. Specifically, the radiative decay Ki was 
considered and was shown to be strong relative to the non-strange contri- 
butions only in a limited kinematic domairP^. 

A final cautionary remark is in order here. When the rate spectra are 
studied at photon energies above 1 GeV, one must keep in mind that 
hadronic form factors have not been implemented in most of the rate cal- 
culations up to now (an exception is the Kapusta, Lichard and Seibert 
calculationiii^ which estimated the form-factor effect on the tttt chan- 
nel), and could result in a suppression of a factor of 2 or more at higher 
photon energies. This is where the exchanged meson goes further off shell 
and brings forward possibly large form-factor effects. Advances in this di- 
rection will be important. 

2.3.4. Medium Effects 

In terms of higher-order effects, corrections to these rates come in at least 
two forms. First, there are off-shell effects which can be conceptualised by 
dressing the propagators and vertices for the internal hadron species in the 
general photon self-energy structure. These are the so-called form factors 
mentioned earlier in cautionary remark. Second, and beyond this, there are 
bona fide medium effects (finite temperature and density effects, e.g. width 
smearing and pole mass adjustments) that could be quite important. The 
typical pursuit in studies of medium modifications is to investigate the ef- 
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fects of dramatic collision broadened vector meson spectral distribution^^ 

and/or the dropping of the rho mass according to the so-called Brown- Rho 

de- 
scaling or some other ansatz. Several a uthor s have studied various pieces 

of the overall medium dependenceJ^(Z32SE2 The trends are the follow- 
ing. While the in-medium vector meson widths are expected to be rather 
large, the effect on photon production is not too significant. This makes 
sense since the vector spectral distributions contribute to photon produc- 
tion only as an integral over the specific distribution — and smearing the 
distribution does not affect the normalisation. Mass shifts, on the other 
hand, have been shown to affect t he rates by anywhere from a factor of 3 
up to an order of mag The results are too model dependent to 

make specific concluding statements at present. 



2.3.5. Alternative Approach: Chiral Reduction Formulae 

Instead of computing photon production rates using a channel-by-channel 
assessment, Steele, Yamagishi, and Zahed used chiral reduction formulae 
together with a virial expansion and they came forward with photon and 
dilepton emission rate estimates. The general idea is that the invariant pro- 
duction rate is proportional to the trace over a complete set of hadronic 
states of the hadronic (Boltzmann weighted) Hamiltonian convoluted with 
a current-current correlator. The hadronic part of the correlator is written 
as a virial type expansion truncated in a particular way. The expansion 
coefficients are constrained by various general arguments, e.g. broken chi- 
ral symmetry, unitarity, and gauge invariance and also, when available, 
constrained by observed spectra: electroproduction, r decay, radiative pion 
decay, and so on. The thermal photon emission estimates in this approach 
tend to be l arger t han those using an effective Lagrangian approach by a 
factor of 2- j' '^ l ^'^ l At present, this might be the honest theoretical error 
bar in the rate estimates even after a decade of model calculations. Progress 
continues especially with effective theories in the hadronic matter converg- 
ing with results from models firmly rooted in QCD as the fundamental 
degrees of freedom as discussed in the next section. This is the so-called 
duality of hadronic matter and quark matter at the phase boundary that 
one expects. 

At this stage in the discussion it is somewhat premature to integrate 
the photon production rates over a space-time evolution, which would then 
facilitate a comparison with experiment, because radiation from partonic 
matter has not yet been discussed. So, before considering photon yields 
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from nuclear collisions and making contact with data, the partonic contri- 
butions to electromagnetic radiation will be presented, and then yields will 
be discussed. 

3. Radiation from Partons 

It is of great theoretical importance to establish the production rates of 
electromagnetic radiation for many-body systems beyond the deconfine- 
ment phase boundary of nuclear matter. A model is employed whereby the 
matter is assumed to be fully in the partonic phase. Whilst experimen- 
tal verification of an unequivocal identification of thermalised quark-gluon 
plasma is still forthcoming, it is the appropriate picture with which to 
work as a baseline. The general formalism established for photon produc- 
tion rates from hadronic matter in Sect. El is generic to all quantum field 
theories and is thus equally valid for partonic degrees of freedom. And yet, 
the massless nature of the up and down quarks requires special attention. 
Calculational tools known as hard-thermal-loop (HTL) methods have been 
applied to handle infrared singularities. Independent of the experimental 
advancements then, it would already be important to establish quark mat- 
ter radiative emissivities. Since heavy-ion experiments at the CERN SPS 
and at RHIC have most likely probed into small areas of the deconfined 
region in the nuclear matter phase diagram, there is further motivation, 
and indeed some urgency, for theoretical investigations to converge and to 
report emission rates. Therefore, the status of theory for photon produc- 
tion from finite temperature quark matter is discussed below and a separate 
section is devoted to dilepton production. 

3.1. Photons 

The imaginary parts of one-loop contributions to the photon self-energy, 
obtained with appropriate cuts, are identically zero due to vanishing phase 
space. Certain two-loop diagrams give nontrivial contributions. Cutting 
rules provide a bridge between kinetic theory and field theory where in 
fact, a mapping has been established^. The result of cutting two-loop dia- 
grams gives QCD processes of the types qq and qg qj or qg q^. 
These processes, as well as bremsstrahlung pr ocesses w ere studied using 
perturbative matrix elements two decades agd ^ 1 ^'^ * ^ I The results were 
unfortunately infrared unstable {i.e. the rates diverged as the quark mass 
tended to zero). Significant improvement came when t he "annihilation" and 
"Compton" processes were analysed by Kapusta et and Baier et 
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using resummation techniques of Braaten and 

PisarsklS3H3. 

The basic idea 

behind the resummation technique, or the so-caUed hard-thermal-loop ap- 
proximation, is that weak coupling at high temperatures allows a separation 
of scales, and a separation of the rate into soft (quark momentum ^ gT or 
smaller) and hard (quark momentum T or larger) contributions. The soft 
contribution can be computed with an appropriately dressed quark prop- 
agator in the one-loop photon self-energy, while the hard contribution can 
be computed using perturbative methods and kinetic theory. In each result, 
the separation scale appears as a sort of regulator. When the soft plus hard 
contributions are collected together and added, the result is independent of 
the separation scale, and of course also independent of quark mass since it 
was set to zero from the beginning. The exact result can be established only 
numerically. However, using an approximation which is valid for E^/T ^ 
1, a simple pocket formula has been proposed. At the time, this result was 
thought to be complete to order aas ■ Specialising to two quark flavours the 
result is 

A value of as = 0.4 {g^ = 5) is used and a "1" is added to the argument 
of the logarithm when plotting as suggested by Kapusta et al. to more 
closely match the exact numerical result for photon energies of the order 
of the temperature. This also ensures the rate is always positive. These 
results were subsequently generalised to finite quark chemical potential and 
also applied to chemical non-equilibrium systems. For a discussion of these 
effects see [HH| and references therein. 

Having photon emission rates from QCD free from infrared instability 
ailments represented significant advancement and was at the time, thought 
to be the complete lowest order result. After all, the two loop contribu- 
tions to the self-energy (that is, the dressed two-loop contributions, which 
actually contain arbitrarily many loops) seem naively to contribute to pho- 
ton production at 0{aa1). They specifically correspond to bremsstrahlung 
processes and annihilation with scattering. The extra vertices would intro- 
duce an extra power of g^ as compared with the one-loop result. However, 
Aurenche et al. showed that the two-loop HTL contributi on i s curiously 
not of higher order, but instead contributes to order aas tocP. This owes 
essentially to a coUinear singularity when the exchanged gluon is soft. The 
resummed gluon propagator introduces a g^ in its denominator which can- 
cels the "extra" g^ from the additional vertices. The overall contribution 
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to the HTL for this category of two-loop diagrams is the same (lowest) 
order in aa^. In terms of the kinetic theory analog, these correspond to 
such processes as qq gcpi, gq gqj, and qqq — > (77, or qqg ~* gj (and 
still others with antiquaries). Two- and even three- loop contributions were 
shown to contribute to lowest order, and the rates continued to rise! 

It is fair to say that after these features were pointed out by Aurenche 
et ai, the situation appeared to signal a breakdown in perturbation the- 
ory for finite temperature QCD. However, it has been shown recently by 
Arnold, Moore and Yaffe that as long as ^ gT, there is sufficient cancel- 
lation due to many-b o dy e ffect s so that the lowest-order rate is identifiable 
and fully under control^. This remarkable result represents a significant 
advancement in this field. There were unfinished details within the top- 
ics of bremsstrahlung^^, magnetic mass^-*^ and coherence effect that 
Arnold et al. resolved by analysing multiple-loop ladder diagrams which 
introduce multiple scattering interference effects of Landau-Pomeranchuk- 
Migdal (LPM)^"^ -231 A digression will not be taken here to reproduce the 
lengthy and specialised argument, but the result is the following. The 
suppression is sufficient to regulate the rates because 1-loop, 2-loop and 
multiple-loop diagrams can be consistently resummed to give a finite rate! 
The efforts of many people over a decade of work have produced a com- 
plete photon production calculation from QCD to lowest order 0{aas). 
The simple expression below parametrises the exact numerical solution for 
two quark flavours. 

dR 5 aag ^ 1 
^dV^"9 3^ e^-/^ + 1 
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where 



rem 



C2^2 ^ 0M1{T/E^) ~ 0.3615 + l.Ole-^-^^^^^^ 

0.633 ln(12.28 + {T/E^)) Q.l^A{E^lT) 



,(29) 



Cbren. + annih - (i5^/T)3/2 + ^ I + [E, / 1&.27T) '^^^^ 

Results are shown in Fig.^| where a value for the strong coupling = 
0.4 {g^ = 5) is used as before, and superimposed onto the total hadron rate 
discussed previously. The striking feature is that after the dust has settled 
on the QCD calculations, with HTL to 1-, 2-, and even multiple-loop order, 
with LPM effects carefully included, the QCD rate at fixed temperature is 
once again the same as the hot hadronic gas rate. The QGP and the hadron 
gas seem to "shine just as brightly" . 
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Fig. 18. Photon emission rate for QGP from Kapusta et al. in Ref. (labelled 
Quarks-II) and Arnold et al. in Ref. 1901 (labelled Quarks— I). Quarks-II includes Comp- 
ton and annihilation, while Quarks— I includes in addition, bremsstrahlung and certain 
3— >2 processes. Quarks— I is the "complete lowest order calculation". Temperature is fixed 
at T = 200 MeV. Total hadron contribution is also displayed for comparison purposes. 

3.1.1. Photon Measurements 

Photon experiments using heavy-ion beams are notoriously difficult and 
signals of any kind are already a notable accomplishment. At high energies, 
there are at present two sets of data with which to compare the theory. 
The WA80 collaboration at CERN first reported and discussed their yields 
as absolute measurements, but were later forced to loosen the constraints 
somewhat and suggest upper limits on ly. Their direct photon limits came 
from 200^ GeV + Au collisionJ^SI. Secondly, the WA98 collaboration 
measured direct photons in 158^ GeV ^°®Pb + ^"^Pb collisions also at 
The hope from the onset was to challenge the theory using pro- 
duction rates convoluted with a temperature profile evolving according to 
one of two possible scenarios: 1) the system first comes into equilibrium 
well above the phase boundary and therefore the quark rates contribute 
until such time as the system reaches the mixed phase. Overlapping four 
volumes mean that quarks and hadrons contribute until the latent heat is 
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absorbed fully into a hadronic state, and finally, the hadrons emit until 
freezeout; and scenario 2) where the system reaches a very hot and dense 
hadronic state and simply radiates photons while cooling and eventually 
freezes out. The burning question is which scenario is consistent with the 
measurements? Can either one be ruled out? 

There were several attempts to describe the WA80 results and do just 
that. Shuryak and XionJ^ first used the hadron rates with their version 
(incoherent treatment) of ai meson dynamics included, and their conclusion 
was that the excess photon signal could not be described with a conventional 
expansion scenario. They consequently suggested a long-lived mixed phase 
as a possible explanation. Since the data were later reported as upper limits 
only, the conclusion no longer rested on strong experimental support. 

Srivastava and Sinha applied the quark rates at the 1-loop HTL level 
and the hadron ra tes comparing scenarios (I) with and (II) without a phase 
transition to QGF^SHI. They argued that the data (which later became up- 
per limits) are well described by a scenario where QGP is formed initially. 
Bjorken hydrodynamics was employed with Ti = 203 MeV, Tc — 160 MeV, 
and Tf = 100 MeV for scenario (I) and, Ti = 408 MeV for scenario (II). 
Other models came forward attempting to describe the experimental re- 
sults. For example, Dumitru etoL used a three-fluid hydrodynamics with- 
out and with a phase transitioJ^. They came to similar conclusions, that 
without a phase transition to quark matter, the results were inconsistent 
with experiment. 

Improvements in rate calculations from quark matter brought advance- 
ment also in yield estimates. Two-loop HTL rates were coupled with hy- 
drodynamics, and then later corrected due to numerical errors along the 
wajff^'^. The most recent and corrected comparison of the WA80 upper 
limits to hydrodynamic model estimates are displayed in Fig 1 191 The con- 
clusion is that both scenarios, without and with a phase transition, seem to 
be consistent with the upper limits. A more complete hadronic equation of 
state (EOS) and up-to-date photon rates from quark matter lead to these 
new conclusions. 

With the WA80 results in hand and even anticipating the forthcom- 
ing WA98 data at that time, Cleymans, Redlich and Srivastava-*^^^ used a 
hydrodynamical model, which arguably provided better description of the 
evolution as compared to previous model calculations and, in particular, 
could more completely describe the transverse flow likely to be generated 
at the SPS. The initial QCD rates of Kapusta et al. were used and the 
hadron rates, including the effects of the ai meson, were implemented with 
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Pt (GeV/c) (GeV/c) 

Fig. 19. Upper limits at the 90% confidence level from WA80 o n th e invariant excess 
photon yield per event for the 7.4% a-^^, most central collisions 1^3 as compared with 
hydrodynamical models with (left) and without (right) phase transition 



all hadrons up to 2.5 GeV mass contributing to the equation of state. They 
concluded that while the final yields were not significantly different in a 
QCD plus hadron matter scenario as compared with a fully hadronic pic- 
ture, they argued that the physics seemed to favour the former since in 
the hadronic picture particle densities were beyond anything reasonable for 
hadronic language to be justified. 

Before moving to the WA98 data, one might make the remark that since 
the WA80 results are upper limits rather than measurements, and due to 
the uncertainties in the theoretical production rates and, mostly, with the 
uncertainties in the models for the evolutions of the nuclear systems, no 
definite conclusions can be reached. 

The eagerly anticipated direct photon measurement from ^"^Pb + 
208pb collisions at 158.4 GeV were published in 2000 by the WA98 
collaboratioJ^. The collaboration presented their data as compared with 
several proton-induced reactions at similar energies and scaled up to central 
208pj^ _l_ 208p|-, collisions. For pT > 1-5 GeV, where the signal is strongest, 
there is a clear excess beyond that which is expected from proton-induced 
reactions. In other words, the results are quite suggestive of thermal pho- 
ton emission, or perhaps pre-equilibrium emission. One contribution that is 
non-negotiable in those data is that due to perturbative QCD. It owes its ex- 
istence to collisions during the first instants of the reaction, and should ap- 
pear in pp, pA, and AA measurements. The WA98 data is shown in Fig. 1201 
along with a pQCD estimatJ^. Even though the presence of p QCD ef- 
fects at the energies under discussion here can't be argued againslM^, the 
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pQCD (s"^ = 17.3 GeV, y=0), Wong e 

CTEQ4, <kT>=0.9 (GeV/cf 
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Transverse Momentum (GeV/c) 

Fig. 20. The WA98 real photon measurements as a function of photon transverse mo- 
mentum. The pQCD estimates are from 11021 . 



application to nucleus-nucleus data contains some ambiguities that need 
to be pointed out in order to make progress. Specifically, it is clear that 
some amount of parton intrinsic transverse momentum (neglected in strict 
NLO calculations) s houl d manifest itself. Simple uncertainty principle ar- 
gumen ts su pport thiJ^, and soft gluon emission should increase the value 
furtheM 

However, attempts to extract meaningful values from experi- 
ments have remained inconclusive; for example a recent survey found that 
fixed target data at ISR energies {^/s < 23 GeV) were inconsistenl^^ 
Furthermore, in nucleus-nucleus collisions, a part of the parton transverse 
momentum can be ascribed to multiple soft scattering of the nucleons prior 
to the hard scattering-'^'-'^, and this has to be modeled dynamically and 
independently. It is important to note that, at RHIC, several of those un- 
certainties will be lifted, as measurements of pp, pA, and AA reactions will 
be performed at the same energy with identical detector configurations. 
Bearing all those caveats in mind, a recent studj ^^'^ of E704 and WA98 
data found that (fc^) ~ 1.3 GeV^ could by extracted from pp reactions, 
leaving up to 1 GeV^ for nuclear effects. This analysis is shown in Fig. |^ 
It is clear from this work that photon transverse momenta below 2.5 GeV 
are under-predicted by this pQCD estimate. Also, around this momentum. 
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Fig. 21. The WA98 data with pQCD estimates of transverse parton momentum. From 
IW7I . 



the exact value of the intrinsic transverse momentum ceases to be impor- 
tant. A softer component of the photon spectrum is caUed for, and this wiU 
be discussed shortly. Note that this value for a separation of scale between 
the "hard" and "soft" photon sources also appears i f one fit s the high mo- 
mentum pQCD spectrum to the data, with a K factoi^SIS] it is argued in 
these cited works that the soft component possesses thermal characteristics. 

Srivastava and Sinhe^^ studied mechanisms for excess photon produc- 
tion using an hydrodynamic expansion applied to the ^°*Pb-|- ^°*Pb system. 
The photon emission rate from quarks was input using the result from two- 
loop HTL calculations from Aurenche et a . However, it is probably fair 
to say that those rates have been superseded by the calculations in Ref. [201 , 
which incorporates higher loop topologies and thus LPM effects. Produc- 
tion rates from the hadronic phase were taken from the parametrisation of 
Kapusta et al. plus an incoherent ai-exchange contribution to the process 
np — !■ TTj. The results obtained there are shown in Fig. 1221 The high initial 
temperature in this work is needed to generate a sufficient high transverse 
momentum component of the photon spectrum. In this respect, the WA98 
data has been used to extract phenomenologically an initial radial veloc- 
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kT (GeV) 

Fig. 22. Single photon production in Pb+Pb collisions at relativistic energies. For spe- 
cific model details, consult Ref. 11091 . 



ity profiljnni. 

The result of that study is shown in Fig. Both of those 
theoretical efforts concluded that the excess seemed to be consistent with 
a thermal source of photons at roughly ~ 200 MeV temperature, while de- 
tailed and quantitative conclusions on a partonic scenario versus hadronic 
with strong flow were not definitively reached. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 0.5 I 1.5 2 2.5 3 3.5 4 4.5 

p,(r,eV) p,(r,eV) 

Fig. 23. Total photon yields from quark matter plus hadronic matter (left panel) and 
hadronic matter with medium-modified vector meson properties (right panel). The figure 
is reproduced here from Ref. IllOI . 
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Ruuskanen and collaborators also used hydrodynamics to compare the- 
ory with experiment -'^-'^-'^ ^^"^ ^^'"^ . These workers have challenged hydrody- 
namics to find consistency with not only photon spectra, but also hadron 
and dilepton spectra — all within the same model and simultaneo usly. They 
also insist on reproducing the longitudinal hadron characteristic a I Sev- 
eral equations of state and therefore several expansion scenarios seem to 
describe the photon spectra and hadron spectra equally well. The degen- 
eracy between the different equations of state and initial conditions is not 
lifted empirically, even though the data do require a high density and tem- 
perature initial phase. 




12 3 4 

Pt (GeV) 

Fig. 24. Comparison of the WA98 photon spectrum to the predictions of the UrQMD 
model and the hydrodynamic model at several freezeout temperatures from Ref. |25[. 

A comparison of space-time models in reference to the WA98 data was 
recently carried out by Huovinen, Belkacem, Ellis and Kapustai*^^' wherein 
hydrodynamics and a course-grained UrQMD were used to produce photon 
spectra (dilepton spectra were also calculated within theory and compared 
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to experiment). Notably, the "complete" lowest order photon production 
rate [(^(aas)] from the quark phase was used in this work. The basic con- 
clusions were that UrQMD and hydrodynamics seem to give roughly the 
same qualitative features of expansion and cooling (although quantitatively, 
UrQMD cools more slowly due presumably to viscous and heat conduction 
effects); they give therefore, very similar results for photon production. The 
precise choice of freezeout temperature seemed to be irrelevant, indicating 
that the high temperature part of the evolution dominates photon produc- 
tion. Results are shown here in Fig. |^ The agreement between theory 
and experiment was described by these authors as "excellent", while they 
reminded the reader that the rates have uncertainties and the initial condi- 
tions which were fed into the models are responsible for further uncertainties 
propagating to the final spectra. 

A partial summary of the photon analyses is justified. It is an accurate 
statement that definite conclusions are elusive. Many physical ingredients 
have been invoked in the studies of heavy ion photon data, as seen above 
and in the quoted reviews, but uncertainties in many of those ingredients 
(if not all) preclude a clear interpretation of a signal that relies on a com- 
bination of their effects. But one example is the absence of the chemical 
potentials in hydrodynamics-based approaches. Another is the uncertainty 
in the basic photon rates. However, those uncertainties have narrowed down 
considerably in recent years, and this is true for rates in both the partonic 
and confined sectors. Also, as mentioned previously, the fact of being able 
to access data at the same energy in pp, pA, and AA events will make 
RHIC a fertile testing ground for theoretical models, and should allow the 
community to make more progress in differentiating between them. 

3.2. Dileptons 

The yield of low mass dileptons (M < m^) from thermal quark- antiquark 
annihilation is not expected to be a great competitor of the two-pion an- 
nihilation simply owing to longevity effects in the two phases. The quark 
phase occupies a smaller four volume. Nevertheless, it is useful to assess the 
production rates as a benchmark and then to ask about higher-order cor- 
rections, especially in the medium. For high enough system temperatures, 
and for large enough invariant masses, 99 — *■ 7* — > e~^e~ is considered a 
more significant source and in terms of theory, can be reliably computed in 
a HTL approximation. The lowest order contribution [0(q!°)] in field the- 
ory language corresponds to a one-loop graph with bare quarks occupying 
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the internal lines. The imaginary part of the self-energy describes precisely 
the annihilation process mentioned above. The production rate is roughly 
the square of the density of quarks times the c ross section times the rel- 
ative velocity. These are now textbook formula^^^ so one simply quotes 
the results 



dR 
dlVP 



where the annihilation cross section is 



, , 47r 



(31) 



(32) 



and where Af is an overall degeneracy factor (24 when using two quark 
flavours) and finally, Ki is the modified Bessel function of order 1. Here the 
quark and lepton masses have been set to zero. To compare the resulting 
rate with major hadronic contributors, Fig. 1251 is presented. Already here, 
one sees that the Born term is not negligible and the natural next question 
is the role of higher order contributions. 
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Fig. 25. Thermal production rate for dileptons via lowest order quark-antiquark anni- 
hilation as compared to leading hadronic channels tttt — > e+e~ and KK — > e+e~. The 
temperatures are set to T = 150 MeV (left) and 200 MeV (right). 
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Perturbative corrections to these annihilation rates were considered by 

1 1 Q 

Braaten, Pisarski and Yuan who found that for very soft dileptons at 
rest in the medium (energies ^ 1 GeV) the corrections were orders of mag- 
nitude larger than the Born term. Also, unique structures emerged in these 
corrections owing to Van Hove singularities arising from significant soften- 
ing of the quark dispersion relation in medium. There appears a minimum 
in the medium-modified quark dispersion relation (a plasmino) typically 
at dilepton energies less than that set by the quark mass. While these 
effects are quite intriguing, finite imaginary parts in t he q uark propaga- 
tors and finite three-momentum effects for the dilepton I could dampen 
the peaks i nto u ndetectable artifacts. Also, the softer bremsstrahlung 
contribution J-^^^l might overshine these total annihilation channels. For a 
review of these and other issues for the dilepton channels, see Ref. |24| by 
Rapp and Wambach. 

Quark-antiquark annihilation is of course not the only relevant par- 
ton process for dilepton production. For instance, the 2^2 real photon 
production processes considered previously contribute also to lepton pair 
emission. In addition, there are annihilation processes where one of the in- 
coming partons has already scattered and suffered an off-shell interaction. 
The resulting 3^2 process comes from off-shell annihilation (also called 
annihilation with scattering). Such mechanisms have been shown to dom- 
inate at high enough photon energy. Since this is essentially a many-body 
initial state, formation time considerations and coherence effects for the 
virtual photon suggest once again that mult iple scattering plays an ir npor- 
tant role. Aurenche, Gelis, and Zarakel^^, and together with Moor^^^l 
have applied the HTL technique for lepton pairs with E/T 3> 1 (either 
low mass but high momentum, or high mass) and have shown that the 
two-loop contributions which include bremsstrahlung of a quark and an- 
nihilation with scattering, are free from infrared and coUinear singularity 
effects. When adde d to to t he Born term, the rescattering corrections plus 
the 2^2 processe^^n^^ result in a rate that is somewhat increased as 
compared with just the Born terms. Furthermore, threshold effects at 
— 4m^ (thermal quark mass) are smoothed out, all of which is illustrated 
in Fig. 

The dilepton results from QGP discussed above assume equilibrium and 
and use asymptotic values for such quantities as thermal quark masses and 
screening masses. Strictly speaking, these are only valid at asymptotic val- 
ues of temperature: the assumptions needed for the theoretical machinery 
to remain consistent might actually break down at terrestrial accelerators 
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Fig. 26. The dilepton production rate per unit four momentum at fixed energy as a 
function of photon mass (dilepton mass). Values are fixed at T = 1 GeV, go = 5 GeV, 
as = 0.3 and two quark flavours were considered. The figure is reproduced here from 
Ref. W^ . 



energies. Scenarios more realistic for RHIC and LHC could be studied if al- 
ternative schemes were used to compute masses in nonperturbative circum- 
stances. First steps in the direction o f latt ice evaluations of thermal dilep- 
tons using maximum entropy method^^^ have recently been takenl^. 



4. Predictions 
4.1. Photons 

RHIC has been running, looking primarily at hadronic observables probing 
the later stages of ultrarelativistic nuclear reactions. Electromagnetic spec- 
tra will soon be available which, of course, probes deeper into the fireball 
and indeed cleanly into early stages of the reactions as well. Predictions 
for measurements of electromagnetic signals at RHIC are therefore very 
important. In addition, the Large Hadron Collider (LHC) is only five years 
away! While this number probably needs an appropriate error bar, it will 
soon become crucial to have formulated a set of model estimates for LHC 
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experiments too. A section is devoted here to discussing these sorts of pre- 
dictions. 

As one moves away from SPS systems and energies and goes to RHIC, 
and to LHC energies, there are increasing uncertainties in estimates for the 
initial energy densities. The initial state is very far from under control. But, 
as in all cases, when theory is extrapolated to new territory, the simplest 
estimates are first used to set the scales and subsequent to this, refine- 
ments and various improvements are made. It is in this spirit that photon 
produc tion |yjglds) wer e recently estimated at RHIC and LHC by several 

authorjnniminHini. 

Simple 1-1-1 dimensional models show dominance of 
the QGP over the hadron gas for photon px ^ 3 GeV (RHIC) and roughly 
2 GeV (LHcjT^i. Transverse expansion, which builds up particularly later 
in the hadron phase, makes distinction less clear, but the QGP might still 
outshine the hadron gas. The results for photon production from Ref. jl29| 
are displayed in Fig. 1771 

However, agreement is far from complete on this issue. For instance, 
Hammon et a/. predict that QGP will not be visible at RHIC owing to 
a very strong contribution from prompt photons (a pre-equilibrium source 
which has not been discussed here, and one which is probably not entirely 
under control), while at LHC the situation is different where QGP will 
be visible for limited photon kinematics. They used Tq = 533 MeV (300 
MeV) and To = 0.12 fm/c (0.5 fm/c) for QGP (hadron gas) at RHIC, 
and they used Tq = 880 MeV (650 MeV) and tq = 0.1 fm/c (0.25 fm/c) 
for QGP (hadron gas) at LHC. One could however hope that the prompt 
photons due to pQCD could be measured separately (in pp collis ions at the 
same energy, for example), and subtracted out. Alam et a/! ^^^ l find, using 
less extreme initial condition parameters (lower initial temperatures), that 
thermal photons will be visible for pT < 2 GeV. However, they also find 
that thermal photons from hot hadronic gas populate the high px region 
even fairly strongly. Again, this is due to a strong flow built up later in the 
hadron phase. 

In the face of such lack of agreement, which owes essentially to large 
uncertainties in the initial conditions, in the nature of the expansion, and 
even in the quark and hadron rates themselves, one suggests that it is pre- 
mature to make any definite statement at this point. In other words, the 
theoretical error bar is too large at present to formulate any physics con- 
clusions from photons. And yet on the optimistic side, theory will progress 
when the newest QCD rates and hadronic rates are implemented into a 
dynamical model which attempts to describe the buildup of collective flow 
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Fig. 27. Photon spectra (yields) from SPS, to RHIC, and LHC from Pb+Pb collisions. 
Two cases are shown, an ideal pion gas with =3 (left) and 8 (right) panels. 



in some detail, on a species-by-species basis (viz. heavier species seem to 
flow difl^erently from lighter ones). 

In almost all cases discussed in this work, the dynamical simulations 
used to model the dynamics of nuclear collisions assume some form of 
equilibrium. Many approaches assume both chemical and thermal equi- 
librium, while some only need the latter ingredient. There exists, however, 
a whole class of models that attempt an ah initio rendering of the heavy 
ion reactions. Those are currently the only window one has to the very 
early stages of the collisions, and thus they potentially offer precious in- 
sight on the importance of pre-equilibrium generation of electromagnetic 
radiation. At ultrarelativistic energies, the degrees of freedom that appro- 
priately describe this ph ase a re partonic. A recent prediction of the pho- 
ton yields has been mad cj^^^ l using a version of the parton cascade model 
Along with gauging the importance of the above-mentioned 
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pre-equilibrium effects, this caicuiation involves the apphcation of pertur- 
bative QCD (pQCD) in a domain not necessarily restricted to large mo- 
mentum transfers. The photons there are produced from Compton, annihi- 
lation, and bremsstrahlung processes at the parton level. All lowest-order 
QCD scatterings between massless quarks and gluons are included in this 
model. The obtained photon spectrum for the collision of gold nuclei at 
RHIC is shown in Fig.[2Hl There, the contribution from interactions involv- 
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Fig. 28. Transverse momentum spectrum of photons from central colli sions of gold 
nuclei at RHIC, calculated with the VNI/BMS parton cascade mode^^l (see the text 
for an explanation of the symbols). 



ing at least one primary parton (triangles) is compared with that obtained 
with full binary cascading (diamonds). Most of the photons between 2 and 4 
GeV have their origin in the multiple semi-hard scattering of partons. This 
finding would support the claim that high energy quarks going through a 
quark gluon plasma would yield electromagnetic signatures (see later sec- 
tions). 



4.2. Dileptons 

The plasma signature in the lepton pair channels is expected to manifest 
itself mainly in the so-called intermediate mass sectoJ^(see Section (I2.2|l ^. 
Owing to the large multiplicities germane to the collider conditions, a large 
background will render the extraction of any direct electromagnetic signal 
from the low mass region prohibitively difficult. However, there may be still 
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hope to observe some distortions of the vector meson spectral densities. 
Using a dynamical simulation that accounts for a possible under-saturation 
of the parton chemical abundances, and estimates of the vector self-energies 
in a finite temperature meson gas, the yield on low invariant mass lepton 
pairs was calculated in Ref. |188| and is shown in Fig. 123 It can be seen 
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Fig. 29. Net dilepton yield from an initial plasma phase evolving into a final hadronic 
phase. The full and dash-dotted lines respectively represent cases with and without in- 
medium modifications of the vector meson spectral densities (p, and (fi) 11381 . The 
"cocktail" contribution represents decays from on-shell vector mesons, at the end of the 
hadronic evolution. 



that the in-medium effects translate into a suppression of the p—uj complex, 
and an enhancement below M — 0.65 GeV and above M = 0.85 GeV. The 
broadening of the w is a candidate for experimental ob serva tion. 

Moving to the intermediate mass region, one obtain Jl*^ the results dis- 
played in Fig. 1201 In this calculation, the sensitivity of the results on parton 
equilibrium has been examined, and the reader is invited to consult the rel- 
evant reference for the details. The bottom line, however, is that the quark- 
gluon plasma contribution (as approximated by Born-term qq annihilation) 
does not dip below the Drell-Yan. This being said, the differential invariant 
mass distribution is predicted to be dominated by th e sem ileptonic decays 
of correlated cc pairs, in the intermediate mass regiorPSI. This source has 
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Fig. 30. Intermediate-mass lepton pair spectra at RHIC energies 11381 . Contributions 
from the hadron gas (HG), quark-gluon plasma (QGP), and Drell-Yan (DY) are shown 
individually, along with their sum. 



not been shown in the figures above. However, if the heavy quarks that are 
progenitor for the semileptonic decay loose ene rgy in the strongly interact- 
ing medium, this background will be suppress ed I A direct measurement 
would go a long way in lifting the ambiguitie^^OI 

4.3. Electromagnetic Signatures of Jets 

It is appropriate to discuss the electromagnetic signatures of jets in an 
environment devoted to predictions, as the physics necessary for those to 
exist necessitates high energy and intensity machines such as RHIC and 
the LHC. The fate of high energy jets traversing hot and dense matter is a 
fascinating study, and this whole subfield has become known as that of "jet 
quenching" . The manner in which high energy jets loose energy in a strongly 
inte ractin g medium has been shown to depend on the nature of the medium 
itselPH. Thus, jet tomography is expected to be a sensitive probe of the 
quark-gluon plasma. However, if jets and plasma interact in such a way 
that the jet characteristics are modified, the jet-plasma interactions could 
by the same token lead to the emission of electromagnetic radiation. As 
discussed earlier, the microscopic processes leading to real photon emission 
at the parton level are quark-antiquark annihilation, Compton scattering, as 
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well as bremsstrahlung. Therefore, a fast quark passing through the plasma 
will produce photons by Compton scatte ring with the thermal gluons and 
annihilation with the thermal antiquark a I Those processes are higher- 
order in tts, when compared with photons from initial hard scatterings, but 
they will not form a sub-leading contribution as they correspond to multiple 
scattering (actually, double scattering) which grows with the system size. 

Working out the details, one can show that the rate o f pro duction of 
real photons due to annihilation and Compton scattering ia 



(a)^„^\/■s(■5-4m2) 



X a^^^ {s) ^ ' +{q^q). (33) 



xa(^H.)^^^ + (g-9). (34) 

The fi are parton distribution functions. In order to proceed one may as- 
sume that those may be decomposed as 

f{p) = /thermal (P) + /jct(p) (35) 

where the thermal component is characterised by a temperature T: 
/thermal = exp(— £'/T'). This separation is kincmatically reasonable as the 
jet spectra fall of as a power law and can thus easily be differentiated from 
their thermal counterpart. The phase space distribution for the quark jets 
propagatin g thro ugh the QGP is given by the perturbative QCD result for 
the jet yieldi43J 

/jet P = ^2 .2 \j 36 

gq 7ri?^rp_L d^pj_dy 

xSirj - y)e(T - T,)e(w - T)e(i?i - r) , (37) 

where = 2 x 3 is the spin and colour degeneracy of the quarks, R± is 
the transverse dimension of the system, ~ 1 /p± is the formation time 
for the jet and rj is the space-time rapidity. R{r) is a transverse profile 
function. Tmax is the smaller of the life-time r/ of the QGP and the time Td 
taken by the jet produced at position r to reach the surface of the plasma. 
The boost invariant correlation between the rapidity y and rj is assumec 
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4 6 8 10 12 14 

Pt [GeV] 

Fig. 31. Spectrum dN/d?p^dy of ph oton^ at y = for central collision of gold nuclei 
at \/Smn = 200 GeV at RHIC. Plotteo^^^ is the yield for photons from jets interacting 
with the medium (solid line), direct hard photons (long dashed), bremsstrahlung photons 
(short dashed) and thermal photons (dotted). 



Fig. 1311 contains the resuhs for thermal photons, direct photons due to 
primary processes, bremsstrahlung photons and the photons coming from 
jets passing though the QGP in central collision of gold nuclei at RHIC 
energies. The corresponding results for LHC energies are shown in Fig. 1321 
It is seen that the quark jets passing through the QGP give rise to a large 
yield of high energy photons. This contribution should be absent in pp 
collisions. For RHIC this contribution is the dominant source of photons 
up to p± w 6 GeV. The jet-to-photon conversion falls more rapidly with 
p± than the direct photon yield, similar to a higher twist correction. It is 
clear that this new mechanism for the production of high energy photons 
contributes significantly. In fact, it is the leading source of directly produced 
photons at RHIC in the region pt < 6 GeV/c. 

Very similar considerations should apply to the production of lepton 
pairs, even though the details will of course change. Phase space now allows 
the direct annihilation of a quark and an antiquark into a dilepton. The 
cross section for this process is 

a{M') = ^^,N^{2s + irY.'}' (38) 

/ 

where the sum runs over the flavour of quarks, N^, = 3, and s and Cf stand 
for the spin and the charge of the quark. Using kinetic theory, the reaction 
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Fig. 32. The same as Fig. l31l for central collision of lead nuclei at \/5jvjv = 5.5 TeV at 
LHC. Taken from the same reference. 



rate for the above process can be written as 

R = 



/a(Pa) / — ^/b(pf,)cr(M )Wrel, 



where fi stands for the phase-space distribution of the quark or the anti- 
quark, Pa and pb are their momenta respectively and the relative velocity 
is (for massless quarks) 

EaEb - Pq • Pfc 



W TP 

After some algebra this can be rewritten as 
dR a{AP) 



(40) 



dAP " 2 (27r)6 J ^adXadcj)aXbdxbd(f>bdyadybfa fb 

X 5 [M^ - 2AP iaib cosh{ya - yb) + 2M'^iaXb cos (41) 

where Xa — p?p/M, Xb — p\'/M and ya and yb are the rapidities. The 
integrations over the azimuthal angles yield 

dR _ MV(M2) 



dM 



,f2 



Xa dXa Xb dxb dya dyb fa fb 
1/2 



such that 



4:XaX^ - {2xaXbCosh{ya - yb) ~ ly 



^ ^ 2XaXb COSh(j/a - ^j,) - 1 ^ ^ 



(42) 
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M (GeV) 

Fig. 33. Dilepton spectrunCSl for Pb+Pb at v^sJJi? = 17.4 GeV at SPS. The solid 
line is the result with tq = 0.2 fm/c. The long dashed curves give the results when the 
formation time tq is raised to 0.50 fm/c, thus lowering the temperature. See the quoted 
reference for details. 

< iaM < oo , (43) 
-oo < ya,b < oo . 

When fa and fb are given by a thermal distribution 

/th(p) - cxp{-E/T) =. exp{-pT cosh y/T) , (44) 

the above integral can be performed to obtain Eq. 

Again assuming a Bj or ken-scenario isentropic plasma evolution, one can 
plot results for thermal dileptons, dilcptons from the Drell-Yan process, and 
the dileptons from the passage of quark jets through the plasma for SPS, 
RHIC, and LHC respectively. Those constitute figures 1551 15H and 

At SPS energies, we recover (Fig. 133(1 the well known result that the 
high mass dileptons have their origin predominantly in the Drell-Yan pro- 
cess. Increasing the formation time from 0.20 fm/c to 0.50 fm/c — and 
thus lowering the initial temperature by 100 MeV — drastically alters the 
thermal production (from the dash-dotted curve to the long-dashed one) 
while the yield from the proposed jet-plasma interaction, even though es- 
sentially negligible, is reduced by a factor of w 2 (from the solid line to the 
long-dashed one). 

The jet-plasma interaction starts playing an interesting role at RHIC 
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Fig. 34. Same as Fig.ESJfor central Au+Au at ^/s]^ = 200 GeV at RHIC. 



energies (Fig. I34() , as now the corresponding yield is about only one third 
of the Drell-Yan contribution, and is much larger than the thermal contri- 
bution. Again lowering the initial temperature (now by about 150 MeV) by 
increasing the formation time to 0.50 fm/c further enhances the importance 
of the yield due to jet-plasma interaction. This production is of the same 
order as that attributed to secondary-secondary quark-antiquark annihila- 
tion in a dilepton produ ction calculation done using an earlier version of 
the parton cascade mode^^. 

The much higher initial temperatures likely to be attained at the LHC 
and the much larger (mini)jet production lead to an excess of high mass 
dileptons from jet-plasma interactions which can be an order of magni- 
tude greater (at M = 10 GeV) than that due to the Drell-Yan pro- 
cess. Again, reducing the initial temperature by raising the formation time 
to 0.50 fm/c reduces the jet- plasma yield by about a factor of 2 while 
the thermal yie ld is re duced far more. Recall that at LHC energies sev- 
eral calculation^3n^JL148 j^ave reported a thermal yield larger than the 
Drell-Yan production. It is found that the jet-plasma interaction enhances 
the high mass dilepton production considerably. The calculations out lined 
above can be repeated for a plasma that is not in chemical equilibriuir^^l 
In this case, conclusions similar to the ones reached above are obtained. 
Summarising, it appears that a unique source of high mass dileptons is 
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Fig. 35. Same as Fig.ESJfor central Pb+Pb at ^/s]^ = 5.5 TeV at LHC. 

generated by the passage of quark jets through the quark-gluon plasma. 
The contribution is seen to be the largest at LHC energies, moderate at 
RHIC energies, and negligible at SPS energies. The measurement of this 
radiation could then be added to the list of QGP signatures, as well as 
providing a proof of existence for the conditions suitable for jet-quenching 
to occur. 

Finally, even though this is not a "direct" plasma signal, it is worth 
mentioning that electromagnetic radiation can also serve as a versatile jet- 
tag. This is especially useful in environments where the jet is expected to 
loose energy, or to b e quenched out of existen ce. T his statement holds true 
for real photon^^ll as well as for lepton pair^^DI 



4.4. Squeezing Lepton Pairs out of Broken Symmetries 

We have seen in the text above that the electromagnetic radiation mea- 
sured in nuclear collisions is a precious measure of the in-medium photon 
self-energy. However, in a bath of finite temperature and density, new pos- 
sibilities can manifest themselves. In some sense, the medium allows for 
the existence of correlators that vanish identically in the vacuum. This fact 
opens mixing channels that were previously closed. While several studies 
have sought to investigate the in-medium properties of hadrons, their mix- 
ing with other hadrons has up to now received little attention. An exception 
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is the case oi p — ui EDI This specific mixing may be omitted when deal- 
ing with isospin symmetric nuclear matter, as is done here. Also, one will 
concentrate on vector mesons, as they enjoy a privileged relationship with 
electromagnetic signals. First, one describes an exploratory calculation de- 
signed to highlight an eventual signal. The possibility of p — oq mixing is 
explored, via nucleon-nucleon excitations in strongly interacting systems. 
It is shown that this mixing opens up a new channel for the dilepton pro- 
duction and may thus indu ce a.n add itional peak in the cf) region. A similar 
mixing exists with the a — but its electromagnetic signatures are 

much smaller -'^^'^ ^id) 

For the purposes described above, the interaction Lagrangian can be 
written as 



9p 



2mN 



, (45) 



where VP, (p^, (jiaa, Pi and lo correspond to nucleon, a, ao p, and uj fields, and 
Tf, is a Pauli matrix. The values for the coupling parameters are from [1311 . 
The existence of a preferred rest frame essentially creates a new vacuum 
state with quantum numbers different from those of the true vacuum. An 
immediate consequence of this fact is that one can now define a mixed 
correlator involving scalar and vector current operators: (j^jj). This mixed 
correlator is identically zero in the true vacuum. The polarisation vector 
through which the ao couples to the p via the N — N loop is given by 

U^iqo,\q\) ^ 2iga,gp J ^ Tr [G(fc)r^G(fc + g)] , (46) 
where 2 is an isospin factor and the vertex for p ~ N — N coupling is 

r,^.,-^a,.,^. (47) 

In the above G{k) is the in-medium nucleon propag atojna. For the 
sake of simplicity, one first uses the density-dependent and temperature- 
independent propagator. This approximation will be relaxed later. With 
the evaluation of the trace and after a little algebra, Eq. 021 can be cast in 
a suggestive form: 

n.(<Zo,kl) = ^2, (^2™^-^j| ^ ^4_4(fc.,)2 -(48) 
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TT 



Fig. 36. Feynman diagram for the process tt + tj e+e , which proceeds through p—ao 
mixing. 



This leads immediately to two conclusions. First, q^H^ — 0. Second, only 
two components of the polarisation vector survive after the angular inte- 
gration. This will guarantee that only the longitudinal component of the p 
couples to the scalar meson, while the transverse mode remains unaltered. 
Further note that current conservation implies that out of the two nonzero 
components of 11^, only one is independent. The new mixing channel will 
affect the properties of the mesons in medium, i. e. affect their masses and 
spectral function. Those aspects wil l riqt be discussed at length here, but 
details can be found in the literatur^^^^ZI 

The p — ao mixing opens a new channel in dense nuclear matter, n + r] 
e+e^, which may proceed through A^-A'^ excitations. The Feynman diagram 
for this process is shown in Fig. One may then evaluate cross sections 
for the production of lepton pairs. Evaluating the polarisation in the zero- 
temperature limit, the cross section for this process is 



1 1 2 

[M-^ - ruaj' + miriJM) ^AP - 4ml ' "° ' ' ^^^^ 

where Ho is the zeroth component of the expression in Eq. I|46(l . The nu- 
merical values for the couplings and the calculation details are to be found 
in Ref. jl51| . The cross sections are shown in Fig. 133 for two different val- 
ues of the nuclear density. The familiar process tttt e+e" is also plotted 
to set a scale. A noticeable feature of this plot is that the mixing process 
induces a peak at TOoq = 0.985 GeV. This constitutes a genuine in-medium 
effect which is mostly density-driven: this peak would be completely absent 
in vacuum. Furthermore, the contribution to the cross section at the ao 
mass is comparable in magnitude to that of the tttt channel at its peak. 
Calculations of emission rates where the T = simplification was not made 
also support this clairrPSl Those mixing effects will also affect the bulk 
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Fig. 37. Dilepton spectrum induced by vr + tt ^ e + e and tt + r/ ^ e + e 
considering matter-induced p — ag mixing, (a) p=1.5po (b) p—2.5po 



features of in-medium mesonic behavioun^^. A natural question to ask is 
whether those symmetry-breaking effects in the dilepton spectrum should 
have been observed in any of the past or present experiments. It turns out 
that the required baryonic d ensity is too transient to in fluenc e significantly 
this signal at CERN energie^22 However, the HADES^"^ experiment at 
the GSI, in Darmstadt, Germany, has the resolution and the sensitivity to 
explore the appropriate invariant mass range. Those heavy ion reactions are 
performed at lower energies, with respect to those of CERN, and thus lead 
to higher baryonic densities which persist longer. Evidencing contin ued t he- 
oretical interest at these lower energies, there are new advancement J^* ^!-*- 
in theory pertinent for the dilepton measurements made by the DL^, 
and to be made by HADES. 

There is finally another symmetry-breaking with an electromagnetic sig- 
nature that will be mentioned, even though quantitative evaluations are still 
in their preliminary stage. At zero temperature, and at finite temperature 
and zero charge density, diagrams in QED that contain a fermion loop with 
an odd number of photon vertices are cancelled by an equal and opposite 
contribution coming from the same diagram with fermion lines r unnin g in 
the opposite direction, this is the basic content of Furry's theorerrfr^ (see 

also nHHiiisni). 

In the language of operators we note that these diagrams are encoun- 
tered in the perturbative evaluation of Green's functions with an odd num- 
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ber of gauge field operators: 

(0\Af,iAf,2...Af,2n+l\0). 

In QED we know that CA^^C^^ = — A^ , where C is the charge conju- 
gation operator. In the case of the vacuum |0), we note that C|0) = |0), as 
the vacuum is uncharged. As a result 

{0\A^,A^,...A,,^^,\0) = {0\C-'CA^,C-'CA^,...A^,^^,C-'C\0) 
= {0\A,,A,,...A,,^^,\0)i-ir^+^ 
= -(0|A^,A^,...A^,„^jO) = 0. (50) 

In an equilibrated medium at a temperature T, we not only have the 
expectation of the operator on the ground state but on all possible matter 
states weighted by a Boltzmann factor: 

n 

where (3 = 1/T and fj, is a. chemical potential. We are thus calculating the 
expectation in the grand canonical ensemble. Here, C\n) — e"^| —n), where 
I — n) is a state in the ensemble with the same number of antiparticles as 
there are particles in \n) and vice- versa, li fi = i.e., the ensemble average 
displays zero density then inserting the operator C^^C as before, we get 

{n\A^,A^, . . . ^^,„^Jn)e-''^" 

= -{-n\A^,A^,...A^,^_^,\-n)e-^^- (51) 



The sum over 

all states will contain the mirror term (— njA^^Ap^ . . . A^2^^^-^\—n)e~^^" 
with the same thermal weight, hence 



J2{n\A^.A^, . . . A^,„^, \n)e-^^" = 0, (52) 

n 

(the expectation over states which are excitations of the vacuum |0) will 
again be zero as in Eq. I5UI) and Furry's theorem still holds. However, if 
/i 7^ (=4> unequal number of particles and antiparticles ) then 

= -(-n|A^,A^,...A^,„^J-n)e-^(^"-^Q"), (53) 
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the mirror term this time is {—n\A^^Api^^ . . . ^^jn+i I ~ n)e^^^^"^'^^"\ with 
a different thermal weight, thus 

Y,{n\A,,A,, . . .^^,„^Jn)e-/'(^"-^Q") ^ 0. (54) 

n 

The sum over all medium states which leads to the expectation value 
is no longer zero. This represents the breaking of Furry's theorem by the 
medium. Note that this occurs only in media with non-zero density or 
chemical potential. Summarising, if the medium contains a net charge, 
such that it breaks C spontaneously. Green's functions that vanish iden- 
tically in the vacuum (or in a neutral medium) can survive. Making the 
simplest possible extension to QCD, one may replace two of the photons 
with gluons. This enables processes like gg — s- £'^£~, where the gluon fu- 
sion proceeds through a quark (antiquark) loop. This channel is exciting 
for the following reason: it offers a direct electromagnetic signature of early 
gluon populations. This represents pristine information oi i the state of the 
many-body system. Calculations are technically involvec&^, but results 
are finally forthcoming . 

5. Conclusions 

This ends our survey of the use of electromagnetic signals as probes 
of strongly-interacting relativistic many-body dynamics. The supporting 
framework has been relativistic quantum field theory, generalised to finite 
temperatures and densities, in order to formulate very general computa- 
tional tools for estimating production rates from heated and compressed 
nuclear systems. The focus was on establishing within equilibrium circum- 
stances, the number of radiated electromagnetic quanta per unit volume 
per unit time. For then, one can, and did, take the rate and evolve accord- 
ing to some "best guess" scenarios for the expansion dynamics in heavy ion 
collisions realized at facilities at SPS and RHIC. 

Dilepton experimental circumstances resembling those expected at the 
CERN SPS have been modeled in a variety of ways. Comparisons of theory 
and experiment have been fruitful in terms of suggesting a consistent picture 
of modified vector meson spectral properties. These modification are truly 
collective nuclear effects. Tremendous advancement in our understanding 
of the way in which nuclear matter responds when it is forced near the 
phase boundary between hadrons and quarks has come from these pursuits. 
To mention some specific achievements, one notes that the rho spectral 
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function has been essentially measured at finite energy density and has 
been shown to be significantly modified from its vacuum structure. This is 
a significant achievement. 

The photon studies have been exceptionally fruitful too. Theory has 
advanced from a stage where rate estimates where plagued with infrared 
singularities, to first establishing regulated lowest-order results with new 
computational techniques and tools. The hard-thermal-loop approach is 
useful not only for photon production, but other studies in hot gauge theo- 
ries and for a variety of observables. Next, one witnessed an impressive effort 
to understand the photon self-energy up to the many-loop order and includ- 
ing multiple-scattering effects within the medium. The celebrated result is 
a complete lowest order photon production rate from finite temperature 
QCD that is stable and reliable. That too is a nontrivial accomplishment. 

When the QCD rates and hadronic rates are used to predict photon 
yields and then compared with experiments WA80 and WA98 from CERN, 
there is consistency, if not discriminatory features. One can say that the re- 
sults strongly suggest thermal emission from a fireball at roughly 200 MeV 
temperature. Have we observed the QGP? The investigations reported on 
in this work only contain hints of an answer. As mentioned previously, the 
radiation from the partonic phase is present in the analyses, but does not 
constitute a large portion of the overall signal. Fortunately, this state of 
affairs is directly related to probed temperatures and to the space-time vol- 
ume occupied by the plasma. Both those are expected to increase in the 
current and future generations of collider experiments. The only indirect 
proof is that, in many dynamical simulations, it is unavoidable for the initial 
phase to be elsewhere in the phase diagram than in the deconfined region. 
This is a consistency requirement brought about by our current knowledge 
of the equation of state. This however will change and evolve, especially 
with the experience gained with non-perturbative approaches which will 
in turn guarantee a better focused picture of the quasiparticle nature of 
the partonic sector. Have we observed chiral symmetry restoration? The 
approach to chiral symmetrv is closely related to the properties of the in- 
medium spectral densitie J I As the constraints of the Weinberg sum rules 
are extended to finite temperatures, they require a degeneracy of the vector 
and axial-vector correlators in the symmetric limit. This demand, however, 
can be satisfied in several wayJI^. Thus, a verdict on the status of chiral 
symmetry restoration is being hindered by the difficulty to access the axial- 
vector correlator unambiguously. More theoretical work needs to be done in 
that respect as well, in order to provide a unified calculation with a credible 
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degree of sophistication. It should finaUy be mentioned that RHIC and the 
LHC will have the intensity to make possible Hanbury-Brown-Twiss inter- 
ferometry me asure ments of direct photons. The theory of this observable is 
well-develope and measurements of the correlation functions are ex- 
pected toplace constraints on the space-time extent of the photon-emitting 
source J " I 

If the estimates brought about by dynamical approaches and by analyses 
of hadrons abundances are reliable, we have just grazed the phase boundary 
of the deconfined sector at the SPS. This assertion is not in conflict with 
the evidence obtained from measurements of electromagnetic observables. 
Then, RHIC, and the LHC, should soon signal bold incursions into a new 
territory. 

Acknowledgements 

It is a pleasure to acknowledge helpful comments from, and discussions with, 
P. Aurcnche, J.-e. Alam, A. Bourque, E. Bratkovskaya, A. Dutt-Mazumder, 
S. Gao, F. Cells, P. Jaikumar, B. Kampfer, J. Kapusta, V. Koch, I. Kvas- 
nikova, G. D. Moore, R. Rapp, D. Srivastava, O. Teodorescu, Y. Tserruya, 
and S. Turbide. C. G. thanks J. Bruce, J. Entwistle, J. Pastorius, Ch. 
Squire, and C. Willis for inspiration. This work has been supported in part 
by the Natural Sciences and Engineering Research Council of Canada, in 
part by the the Fonds Nature et Technologies of Quebec, and in part by 
the National Science Foundation under grant number PHY-0098760. 

References 

1. F. Karsch and E. Laermann, this volume; F. Karsch, AIP Conf. Proc. 631, 
112 (2003). 

2. K. Kanaya, Nucl. Phys. A 715, 233c (2003). 

3. Z. Fodor and F. Karsch, JHEP 03, 014 (2002). 

4. M. Stephanov, K. Rajagopal, and E. Shuryak, Phys. Rev. Lett. 81, 4816 
(1998). 

5. See, for example, K. Rajagopal and F. Wilczek, At the Frontiers of Particle 
Phystcs/Handbook of QCD, M. Shifman ed., (World Scientific, Singapore, 
2000), and references therein. 

6. Prashanth Jaikumar, Ralf Rapp, and Ismail Zahed, Phys. Rev. C 65, 055205 
(2002). 

7. E. L. Feinberg, Nuovo Cimento 34A, 391 (1976). 

8. L.D. McLerran and T. Toimela, Phys. Rev. D 31, 545 (1985). 

9. C. M. Ko, L. H. Xia, and P. J. Siemens, Phys. Lett. B231, 16 (1989). 
10. H. A. Weldon, Phys. Rev. D 42, 2384 (1990). 



1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3'v8 



62 Electromagnetic radiation from relativistic nuclear collisions 

11. C. Gale and J. I. Kapusta, Nucl. Phys. B357, 65 (1991). 

12. E. S. Pradkin, Proc. Lebedev Physics Institute 29, 6 (1965). 

13. Sec, for example, Michael E. Peskiri and Daniel V. Schroeder, An Introduction 
to Quantum Field Theory, (Addison- Wesley, Reading, 1995), and references 
therein. 

14. Charles Gale and Peter Lichard, Phys. Rev. C 49, 3338 (1994). 

15. G. Chanfray and P. Schuck, Nucl. Phys. A555, 329 (1993). 

16. G. Chanfray, R. Rapp, and J. Wambach, Phys. Rev. Lett. 76, 368 (1996). 

17. R. Rapp, G. Chanfray, and J. Wambach, Nucl. Phys. A617, 472 (1997). 

18. J. J. Sakurai, Currents and Mesons, (University of Chicago Press, Chicago, 
1969). 

19. F. Klingl, N. Kaiser, and W. Weise, Z. Phys. A 356, 193 (1996). 

20. H. B. O'Comiel, B. C. Pearse, A. W. Thomas, and A. G. Williams, Prog. 
Part. Nucl. Phys. 39, 201 (1997). 

21. Ralf Rapp and Charles Gale, Phys. Rev. C 60, 024903 (1999). 

22. Sangyong Jeon and Paul. J. Ellis, Phys. Rev. D 58, 045013 (1998). 

23. W. Gassing and E. L. Bratkovskaya, Phys. Rept. 308, 65 (1999). 

24. R. Rapp and J. Wambach, Adv. Nucl. Phys. 25, 1 (2000). 

25. P. Huovinen, M. Belkacem, P.J. Ellis, and J.I. Kapusta, Phys. Rev. C 66, 
014903 (2002). 

26. G. Agakichiev et al, Phys. Rev. Lett. 75, 1272 (1995); Phys. Lett. B422, 

405 (1998). 

27. M. Masera et al., Nucl. Phys. A590, 93c (1995). 

28. G. E. Brown and Mannque Rho, Phys. Rept. 269, 333 (1996). 

29. See, for example, C. Gale, Nucl. Phys. A698, 143 (2002), and references 
therein. 

30. Johannes P. Wessels et al. (CERES collaboration), Nucl. Phys A715, 262c 
(2003). 

31. E. Shuryak, Phys. Lett. B79, 135 (1978). 

32. A. L. S. Angelis et al., Eur. Phys. J. C 13, 433 (2000). 

33. M. C. Abreu et al., Phys. Lett. B368, 230 (1996); Nucl. Phys. A610, 331 

(1996). 

34. C. Y. Wong and Z. Q. Wang, Phys. Lett. B367, 50 (1996). 

35. C. Spieles et al., Eur. Phys. J. C 5, 349 (1998). 

36. Z. Lin and X. N. Wang, Phys. Lett. B444, 245 (1998). 

37. L. Capelli, PhD thesis, Universitc Claude Bernard, 2001. 

38. P. Levai, B. Miiller, and X. N. Wang, Phys. Rev. C 51, 3326 (1995). 

39. P. Bordalo et al., Nucl. Phys. A661, 538c (1999). 

40. G. Q. Li and C. Gale, Phys. Rev. Lett. 81, 1572 (1998); Phys. Rev. C 58, 

2914 (1998). 

41. L. Xiong, E. V. Shuryak, and G. E. Brown, Phys. Rev. D 46, 3798 (1992); 
C. Song, C. M. Ko, and C. Gale, Phys. Rev. D 50, R1827 (1994); J. K. Kim 
et al., Phys. Rev. D 53, 4787 (1996). 

42. Particle Data Group, Phys. Rev. D 66, 010001 (2002). 

43. S. Gao and C. Gale, Phys. Rev. C 57, 254 (1998). 

44. L. M. Barkov al., (OLYA collaboration), Nucl. Phys. B256, 365 (1985). 



1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3'v8 



Electromagnetic radiation from relativistic nuclear collisions 63 

45. D. Bisello et al, (DM2 collaboration), Phys. Lett. B220, 321 (1989). 

46. M. E. Biagini, S. Dubnicka, E. Etim, snd P. Kolar, Nuovo Cimento A 104, 

363 (1991). 

47. R. R. Akhmetshin et al, (CMD-2 collaboration), Phys. Lett. B364, 199 
(1995). 

48. D. Bisello et al, (DM2 collaboration), Z. Phys. C 39, 13 (1988). 

49. P. M. Ivanov et al, (OLYA collaboration), Phys. Lett. B107, 297 (1981). 

50. K. Hagliu and C. Gale, Phys. Rev. D 52, 6297 (1995). 

51. S. I. Dolinsky et al, (ND collaboration), Phys. Rept. 202, 99 (1991). 

52. N. Albrecht al, (ARGUS collaboration), Phys. Lett. B185, 223 (1987). 

53. H. Gomm, 6. Kaymakcalan, and J. Schechter, Phys. Rev. D 30, 2345 (1984). 

54. D. Bisello et al, (DM2 Collaboration), Proceedings of Hadrons '91, (College 
Park, 1991). 

55. Z. Huang, Phys. Lett. B361, 131 (1995). 

56. L Kvasnikova, PhD thesis, McGill University, 2001. 

57. Z. Huang, private communication. 

58. J. D. Bjorkcn, Phys. Rev. D 27, 140 (1983); R. C. Hwa and K. Kajantic, ibid. 
32, 1109 (1985): M. Gyulassy and T. Matsui, Phys. Rev. D 29, 419 (1984). 

59. J. Cleymans, K. Redlich, and D. K. Srivastava, Phys. Rev. C 55, 1431 (1997); 
Phys. Lett. B420, 261 (1998). 

60. loulia Kvasnikova, Charles Gale, and Dinesh Srivastava, Phys. Rev. C 65, 
064903 (2002). 

61. K. Gallmeister, B. Kampfer, O.P. Pavlenko, and C. Gale, Nucl. Phys. A688, 
939 (2001). 

62. R. Rapp and E. Shuryak, Phys. Lett. B473, 13 (2000). 

63. D. K. Srivastava, B. Sinha, L Kvasnikova, and C. Gale, Nucl. Phys. A698, 

432c (2002). 

64. M. C. Abreu et al (NA50 Collaboration), Eur. Phys. J. C 14, 443 (2000). 

65. D. K. Srivastava, Phys Rev. C 64, 064901 (2001); J.-Y. OUitrault, Phys. Lett. 
B273, 31 (1991). 

66. J.I. Kapusta, P. Lichard, and D. Seibert, Phys. Rev. D 44 2774 (1991); 
Erratum-ibid. D47, 4171 (1993). 

67. P. Singer, Phys. Rev. 130 2441 (1963); 161 1694 (1967). 

68. H. Nadeau, J.L Kapusta, and P. Lichard, Phys. Rev. C 45, 3034 (1992); 
Erratum Phys. Rev. C 47, 2426 (1993). 

69. L. Xiong, E. Shuryak and G.E. Brown, Phys. Rev. D 46 3798 (1992). 

70. C. Song, Phys. Rev. C 47 2861 (1993). 

71. Ralf Rapp, Simon Turbide, and Charles Gale, in preparation. 

72. M. Halasz, J.V. Steele, G. Li, G.E. Brown, Phys. Rev. C 58 365 (1998). 

73. K. Haglin, Phys. Rev. C 50 1688 (1994). 

74. K.L. Haglin, Nucl. Phys. A 584 719 (1995). 

75. G.E. Brown and M. Rho, Phys. Rev. Lett. 66 2720 (1991). 

76. S. Sarkar, J. Alam, P. Roy, A.K. Dutt-Mazumder, B. Dutta-Roy, and B. 
Sinha, Nucl. Phys. A 634 206 (1998). 

77. C. Song and G. Fai, Phys. Rev. C 58 1689 (1998). 

78. P. Roy, S. Sarkar, J. Alam, and B. Sinha, Nucl. Phys. A 653 277 (1999). 



1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3'v8 



64 Electromagnetic radiation from relativistic nuclear collisions 

79. J.V. Steele, H. Yamagishi, and I. Zahed, Phys. Lett. B 384 255 (1996). 

80. J.V. Steele, H. Yamagishi, and I. Zahed, Phys. Rev. D 56 5605 (1997). 

81. H.A. Woldon, Phys. Rev. D 28, 2007 (1983); R. L. Kobos and G. W. Scmcnoff, 
Nucl. Phys. B260, 714 (1985); S. M. H. Wong, Phys. Rev. D 64, 025007 
(2001); A. Majumder and C. Gale, Phys. Rev. C 65, 055203 (2002). 

82. E.V. Shuryak, Yad. Fiz. 28, 796 (1978). 

83. K. Kajantie and H.I. Miettinen, Z. Phys. C 9, 341 (1981). 

84. R.C. Hwa and K. Kajantie, Phys. Rev. D 32, 1109 (1985). 

85. R. Baler, H. Nakkagawa, A. Niegawa, and K. Redhch, Z. Phys. C 53, 433 
(1992). 

86. R.D. Pisarski, Nucl. Phys. 309, 476 (1988); Phys. Rev. Lett. 63, 1129 (1989). 

87. E. Braaten and R.D. Pisarski, Nucl. Phys. B 337, 569 (1990). 

88. T. Pcitzmann and M. Thoma, Phys. Rcpt. 364, 175 (2002). 

89. P. Aurenche, F. Gehs, H. Zaraket, and R. Kobes, Phys. Rev. D 58, 085003 
(1998). 

90. P. Arnold, G.D. Moore, and L.G. Yaffe, JHEP 0112, 009 (2001). 

91. P. Aurenche, F. Gclis, H. Zaraket, Phys. Rev. D 61, 116001 (2000). 

92. P. Aurenche, F. Gelis, H. Zaraket, Phys. Rev. D 62, 096012 (2000b). 

93. L.D. Landau and L Pomeranchuk, Dokl. Akad. Nauk Ser. Fiz. 92, 535 (1953). 

94. A.B. Migdal, Phys. Rev. 103, 1811 (1956). 

95. R. Albrecht et al.. The WA80 Collaboration, Phys. Rev. Lett. 76, 3506 

(1996). 

96. M.M. Aggarwal et al.. The WA98 Collaboration, Phys. Rev. Lett. 85, 3595 
(2000). 

97. E.V. Shuryak and L. Xiong, Phys. Lett. B 333, 316 (1994); Erratum, ibid., 
B 339, 334 (1994). 

98. D.K. Srivastava and B. Sinha, Phys. Rev. Lett., 73, 2421 (1994). 

99. A. Dumitru, U. Katscher, J. A. Maruhn, H. Stocker, W. Greiner, and D.H. 
Rischke, Phys. Rev. C 51, 2166 (1995). 

100. D.K. Srivastava and B. Sinha, Eur. Phys. J. C 12 109 (2000); Erratum, ibid, 
C 20 397 (2000). 

101. J. Cleymans, K. Redlich, and D.K. Srivastava, Phys. Rev. C 55, 1431 (1997). 

102. C. Y. Wong and H. Wang, Phys. Rev. C 58, 376 (1998). 

103. P. Aurenche, R. Baier, M. Fontannaz, and D. SchifF, Nucl. Phys. B286, 509 
(1987); Nucl. Phys. B297, 661 (1988). 

104. J. C. CoUins, D. E. Soper, and G. Sterman, Nucl. Phys. B250, 199 (1985). 

105. P. Aurenche et al, Eur. Phys. J. C 9, 107 (1999). 

106. Geoffrey T. Bodwin, Stanley J. Brodsky, and G. Peter Lepage, Phys. Rev. 
D 39. 3287 (1989). 

107. A. Dumitru, L. Frakfurt, L. Gerland, H. Stocker, and M. Strikman, Phys. 
Rev. C 64, 054909 (2001). 

108. K. Gallmeister, B. Kampfer, and O. P. Pavlenko, Phys. Rev. C 62, 057901 
(2000). 

109. D.K. Srivastava and B. Sinha, Phys. Rev. C 64 034902 (2001). 

110. J.-e Alam, S. Sarkar, T. Hatsuda, T.K. Nayak, and B. Sinha, Phys. Rev. C 
63, 021901R (2001). 



1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3v8 



Electromagnetic radiation from relativistic nuclear collisions 65 



111. J. SoUfrank, P. Huovinen, M. Kataja, P.V. Ruuskanen, M. Prakash, and R. 
Venugopalan, Phys. Rev. C55, 392 (1997). 

112. P. Huovinen, P.V. Ruuskanen, and J. SoUfrank, Nucl. Phys. A 650, 227 
(1999). 

113. P. Huovinen, P. V. Ruuskanen, and S. S. Rasanen, Phys. Lett. B 535, 109 
(2002). 

114. C. Y. Wong, Introduction to Heavy Ion Collisions, (World Scientific, Singa- 
pore, 1994). 

115. StefTen A. Bass, Berndt MiiUer, and Dinesh K. Srivastava, Phys. Rev. Lett. 
90, 082301 (2003). 

116. K. Geiger and B. MuUer, NucL Phys. B369, 600 (1992). 

117. S. A. Bass, B. MiiUer, and D. K. Srivastava, Phys. Lett. B 551, 277 (2003). 

118. E. Braaten, R.D. Pisarski, and T.C. Yuan, Phys. Rev. Lett. 64, 2242 (1990). 

119. S. M. H. Wong, Z. Phys. C 53, 465 (1992). 

120. K. Haghn, C. Gale, and V. Emel'yanov, Phys. Rev. D 47, 973 (1993). 

121. P. Aurenche, F. Gehs, and H. Zaraket, JHEP 0207 063 (2002). 

122. P. Aurenche, F. Gehs, G.D. Moore and H. Zaraket, JHEP 0212 006 (2002). 

123. T. Ahherr, P.V. Ruuskanen, NucL Phys. 380, 377 (1992). 

124. M. H. Thoma and C.T. Traxler, Phys. Rev. D 56, 198 (1997). 

125. For a review of this technology, see M. Jarrel and J.E. Gubernatis, Phys. 
Rept. 269, 133 (1996). 

126. I. Wetzorke, hep-lat/0305012 M. Asakawa, T. Hatsuda, and Y. Nakahara, 
Nucl. Phys. A 715, 863c (2003). 

127. N. Hammon et al, Phys. Rev. C 57, 3292 (1998). 

128. D.Y. Peressounko and Y.E. Pokrovsky, Nucl. Phys. A 669, 196 (2000). 

129. F.D. Steffen and M. H. Thoma, Phys. Lett. B 510, 98 (2001) 

130. A. K. Dutt-Mazumder, B. Dutta-Roy, and A. Kundu, Phys. Lett. B399, 
196 (1997); W. Broniowski and W. Florkowski, Phys. Lett. B440, 7 (1998). 

131. R. Machleidt, Adv. Nucl. Phys. 19, 198 (1989). 

132. B. Serot and J. D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 

133. S. A. Chin, Ann. Phys. (N.Y.), 108, 301 (1977). 

134. K. Saito, K. Tsushima, A. W. Thomas, A. G. Wilhams, Phys. Lett. B433, 
243 (1998). 

135. H. A. Weldon, Phys. Lett. B274, 133 (1992); G. Wolf, B. Friman, and M. 
Soyeur, Nucl. Phys. A640, 129 (1998). 

136. O. Teodorescu, A. K. Dutt-Mazumder, and Charles Gale, Phys. Rev. C 63, 
034903 (2001). 

137. O. Teodorescu, A. K. Dutt-Mazumder, and Charles Gale, Phys. Rev. C 66, 
015209 (2002). 

138. RaU Rapp, Phys. Rev. C 63, 054907 (2001). 

139. E. Shuryak, Phys. Rev. C 55, 961 (1997); Z. Lin, R. Vogt, and X.-N. Wang, 
ibid 57, 899 (1998); B. Kampfer and O. P. Pavlenko, Phys. Lett. B 391, 185 
(1997). 

140. http://na60.web.cern.ch/NA60/ 

141. Miklos Gyulassy, Ivan Vitev, Xin-Nian Wang, and Ben-Wei Zhang, this 
volume. 



1, 2008 21:17 



WSPC/Trim Size: 9in x 6in for Review Volume 



qgp3v8 



66 Electromagnetic radiation from relativistic nuclear collisions 



142. Rainer J. Fries, Berndt Miiller, and Dinesh K. Srivastava, Phys. Rev. Lett. 
90, 132301 (2003). 

143. Z. Lin and M. Gyulassy, Phys. Rev C 51 217 (1995). 

144. Dinesh K. Srivastava, Charles Gale, and Rainer J. Fries, Phys. Rev. C 67, 
034903 (2003). 

145. K. Geiger and J. L Kapusta, Phys. Rev. Lett. 70, 1920 (1993). 

146. J. Kapusta, L. D. McLerran and D. K. Srivastava, Phys. Lett. B 283, 145 
(1992). 

147. D. K. Srivastava, M. G. Mustafa and B. Miiller, Phys. Rev. C 56, 1064 
(1997). 

148. R. Vogt, B. V. Jacak, P. L. McGaughey and P. V. Ruuskanen, Phys. Rev. D 
49, 3345 (1994); S. Gavin, P. L. McGaughey, P. V. Ruuskanen and R. Vogt, 
Phys. Rev. C 54, 2606 (1996); D. Fein, Z. P. Valerio and L Sarcevic, Phys. 
Rev. C 56, 1637 (1997). 

149. X. N. Wang, Z. Huang, and L Sarcevic, Phys. Rev. Lett. 77, 231 (1996); X. 
N. Wang and Z. Huang, Phys. Rev. C 55, 3047 (1997). 

150. Dinesh K. Srivastava, Charles Gale, and T. C. Awes, Phys. Rev. C, in press. 

151. O. Teodorescu, A. K. Dutt-Mazumder, and C. Gale, Phys. Rev. C 61, 
051901(R) (2000). 

152. O. Teodorescu, PhD thesis, McGiU University 2002. 

153. J. Friese et al, Prog. Part. Nucl. Phys. 42, 235 (1999). 

154. E. L. Bratkovskaya, W. Gassing, M. Effenberger, and U. Mosel, Nucl. Phys. 
A 653, 301 (1999); E. L. Bratkovskaya, W. Gassing, and U. Mosel, ihid. 686, 
568 (2001); E. L. Bratkovskaya, ibid. 696, 309 (2001). 

155. R. J. Porter et al, Phys. Rev. Lett. 79, 1229 (1997); W. K. Wilson et al, 
Phys. Rev. C 57, 1865 (1998). 

156. R. Shyam and U. Mosel, Ihep-ph/OSOSOSS] K. Schekhter et al, 
^ucl-th/0305015 

157. W. H. Furry, Phys. Rev. 51 135 (1937). 

158. G. Itzykson and J. B. Zuber, Quantum Field Theory, McGraw Hill, New 
York, (1980). 

159. S. Weinberg, The Quantum Theory of Fields, Vol. 1, Cambridge University 
Press, (1995). 

160. A. Majumder, PhD thesis, McGiU University 2002; A. Bourque, MSc thesis, 
McGill University 2002. 

161. A. Bourque, A. Majumder, and C. Gale, in preparation. 

162. J. L Kapusta and Edward V. Shuryak, Phys. Rev. D 49, 4694 (1994). 

163. R. Hanbury-Brown and R. Q. Twiss, Nature 178, 1046 (1956); G. Gold- 
haber, S. Goldhaber, W. Lee, and A. Pais, Phys. Rev. 120, 300 (1960); U. 
Heinz and B. Jacak, Ann. Rev. Nucl. Part. Sci. 49, 529 (1999); R. M. Weiner, 
Phys. Rept. 327, 249 (2000). 

164. D. Neuhauser, Phys. Lett. B 182, 289 (1986); D. K. Srivastava and J. L 
Kapusta, Phys. Lett. B 307, 1 (1993); D. K. Srivastava and C. Gale, Phys. 
Lett. B 319, 407 (1993); D. Peressounko, Phys. Rev. C 67, 014905 (2003); 
Jan-e Alam et al, Phys. Rev. C 67, 054902 (2003). 



